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Abstract 


Collisions are one of the key processes shaping planetary systems. Aster- 
oid families are outcomes of such collision still identifiable across our 
solar system. The families provide a unique view of catastrophic dis- 
ruption phenomena and have been in the focus of planetary scientists 
for more than a century. Most of them are located in the main belt, a 
ring of asteroids between Mars and Jupiter. Here we review the basic 
properties of the families, discuss some recent advances, and anticipate 
future challenges. This review pays more attention to dynamic aspects 
such as family identification, age determination, and long-term evolution. 
The text, however, goes beyond that. Especially, we cover the details 
of young families that see the major advances in the last years, and we 
anticipate it will develop even faster in the future. We also discuss the rel- 
evance of asteroid families for water-ice content in the asteroid belt and 
our current knowledge on links between families and main-belt comets. 
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1 Introduction 


It was more than a century ago when Japanese astronomer Kiyotsugu 
Hirayama noticed that some asteroids orbit around the Sun along unusually 
close trajectories (Hirayama, 1918). Hirayama called them groups of asteroids 
probably of common origin, probably not being fully aware of how important 
topic in asteroid science he initiated. At this point, a fruitful area of asteroid 
families related studies has begun. 

Why do scientists care about families? Collisions among small solar system 
bodies are of low probability and therefore unlikely to occur over one's lifetime. 
However, over the age of the system, many such collisions occurred and played 
key roles in shaping our planetary system. 

The dating of asteroid families provides us with a collisional history of the 
main belt (Bottke et al., 2015; Spoto et al., 2015). It provides a way to rewind 
a clock and understand what happened at the beginning of the solar system 
formation and its subsequent evolution. In addition to that, families provide 
a window into the interiors of their parent bodies, otherwise not accessible, 
allowing us to understand better the internal composition and structure of 
the asteroids (Masiero et al., 2015). This, in turn, helps us also study the 
heterogeneity in the protoplanetary disks and mechanisms of planetesimals 
formation. 

Asteroid families also provide unique insights into the forces shaping our 
planetary system. Families provide us with observable evidence of large-scale 
catastrophic impacts, giving us the tools to test impact physics on planetary 
scales (Michel et al., 2015). Furthermore, families have been critical in revealing 
the fingerprints left by mean motion and secular orbital resonances, as well as 
the non-gravitational Yarkovsky effect (Nesvorny et al., 2015; Carruba et al., 
2018a). 

Asteroid collisions create thousands of smaller rocks of different sizes. Since 
their formation, the family members have been subject to different evolution- 
ary processes, including collisional and dynamical evolution. At least some 
well-placed families are likely to supply a number of fragments to the near- 
Earth region, connecting, therefore, these families to changes in the impact 
flux of inner solar system bodies (Heck et al., 2017). 

Finally, the parent bodies of some families were likely water-bearing aster- 
oids. Following the evolution of the fragments formed by collisional disruption 
of such bodies allows estimating a present-day water's distribution and its 
content in asteroids (Hsieh et al., 2018, 2020), relevant in different aspects, 
including the water delivery to the Earth. Asteroid families are natural 
laboratories to study all these phenomena. 

We aim here to put together the results of recent studies on asteroid fam- 
ilies, which help navigate current achievements, open problems, and future 
opportunities. The article focuses on recent advances, achieved mainly since 
the year 2015. We have tried to avoid duplication with similar attempts, espe- 
cially two Asteroids IV chapters on asteroids. Therefore, this review should be 
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considered primarily complementary to these works. However, some overlap- 
ping is unavoidable, mostly when recalling the previous results is needed to 
understand the new works. 

'The paper is organised as follows. In Section 2 we discussed the meth- 
ods for family identification and attribution of new members to the known 
families. Section 3 review the efforts on modelling the long-term evolution of 
asteroid families and related processes. At the same time, Section 4 presents 
the results on special family classes such as young and very young families 
and water-bearing families. Each section ends by discussing future challenges 
and opportunities to study the related open problems on asteroid families. In 
Section 5 we briefly present our conclusions. 


2 Identification 


'This section reviews the existing methods for classifying asteroids into fam- 
ilies and discusses the main challenges that need to be addressed in the 
future. Methods especially suitable for young families, such as the backward 
integration method, are discussed in Section 4.1. 


2.1 Traditional Hierarchical Clustering Method 


A standard approach to identifying asteroid families is the so-called Hierar- 
chical Clustering Method (HCM), introduced by Zappalà et al. (1990). This 
the most often used methodology is extensively described in the literature. In 
particular, the advantages and disadvantages of the HCM are discussed, for 
instance, by Bendjoya et al. (1993) or Nesvorny et al. (2015). Here we only 
briefly review the main idea and properties of the method. 

Detection of asteroid families and identification of their members is typi- 
cally performed in the 3-D space of proper orbital elements: semi-major axis ap, 
eccentricity ep, and inclination ij (Milani et al., 2014; Nesvorny et al., 2015), 
though in principle this could be done also in the space of proper frequencies 
(Carruba and Michtchenko, 2007, 2009). 

'The distances among the asteroids in the proper elements space are usually 
determined using the following metrics proposed by Zappalà et al. (1990): 


5 f Aap V? 
d = nap 1 Ka + 2(Ae,)? + 2(A sinip)?, (1) 


where nap is the heliocentric orbital velocity of an asteroid on a circular orbit 
having the semi-major axis ap = (@p,2 + ap,1)/2, Aap = Gei — Mp2, Aep = 
€p,1 — €p,2 and A sin ij = sin ip 1 — sin ip,2, where the indexes (1) and (2) denote 
the two bodies whose mutual distance is calculated. The distance d is usually 
expressed in meters per second. 

The HCM used to identify asteroid families is based on a single-linkage 
clustering. It works in an agglomerative (bottom-up) fashion so that at each 
step, two clusters that contain the closest pair of elements, not yet belonging 
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to the same cluster, are combined. The main advantage of the single-linkage 
HCM is that it does not preclude apriori any shape of the family. Therefore, it 
allows the detection of irregular clusters, i.e. those which have irregular shapes 
with narrow bridging regions in multidimensional space. However, the main 
advantage is closely related to the main issue with the HCM, the so-called 
chaining effect. That is, the first concentrations naturally tend to incorporate 
nearby groups, thus forming a chain. In the case of nearby families, the chaining 
could prevent separating the families. The chaining effect is because only a 
single pair of points is needed to be close enough to merge two clusters. 


2.2 Removal of interlopers 


The HCM is a statistical technique that, along with real family members, 
unavoidable associates also some interlopers with a family (Migliorini et al., 
1995). The presence of these interlopers complicates any family-related study, 
and their removal is of great importance. Asteroid families are known to 
be mostly homogeneous in the composition (e.g. Parker et al., 2008), and 
therefore the members of a family share similar surface reflectance characteris- 
tics, including the spectra, albedos, and colours. If available, this information 
could be used to discriminate between actual family members and nearby 
background asteroids.! 

The surface reflectance data could be used in different ways. One possible 
method to exploit these data is to apply the HCM in extended space, i.e., 
in the space that includes surface reflectance properties along with the three 
proper elements. Following pioneering work by Nesvorny et al. (2005a), several 
authors employed this strategy. Parker et al. (2008) applied the HCM in four- 
dimensional space, using a linear combination of the Sloan Digital Sky Survey 
(SDSS) colours as the fourth dimension. Carruba et al. (2013) extend this 
approach, adding albedos as the fifth dimension, significantly reducing the 
percentage of known interlopers with respect to other methods. 

Another possible strategy is to separate the main-belt asteroids into two 
populations (typically representing C- and S-type objects) according to their 
colour or albedo values. The HCM is then run to each of these populations 
separately. Masiero et al. (2013a) utilize this technique onto the WISE albedos. 
They were able to link almost 40,000 asteroids to 76 dynamical families and 
identify several new families. 

A more recent example of this technique was applied by Novakovié et al. 
(2017). These authors used available data to separate asteroids from the Pho- 
caea region into C- and S-complex and then applied the HCM method to the 
sub-population of C-type objects. The approach allowed identification of the 
Tamara family, a group of dark asteroids embedded in the dominant population 
of bright and rocky S-type asteroids. 


lAs a matter of completeness, let us mention that some interlopers could also be identified 
based on their position concerning the V-shape of the family (Vokrouhlicky et al., 2006b; Spoto 
et al., 2015). Such objects are located outside the family V-shape lines. Therefore, they are too 
far from the centre of the family to be transported to their present location by the Yarkovsky 
effect. This is illustrated in Section 3.2 for the Massalia family. 
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Both the above-described strategies are very useful but have limitations in 
that they can only be applied to a reduced set of main-belt asteroids for which 
the colours and albedos are obtained. Despite the significant increase of avail- 
able physical data in recent years, the number of asteroids for which these data 
are at our disposal is still several times smaller than the number of objects for 
which proper elements have been computed. As Milani et al. (2014) explained, 
the proper elements contain more comprehensive information than the physical 
data because the latter are available either for significantly smaller cata- 
logues or with lower relative accuracy. For this reason, some authors adopted 
a bit different approach, so that the so-called dynamical families are first 
obtained in the space of proper elements, and available physical data is used 
only posterior to identify interlopers among family members (e.g. Novakovié 
et al., 2011; Milani et al., 2014). While it allows identifying much more fam- 
ily members with respect to the previous approach, the pitfall of this strategy 
is that reflectance characteristics are of limited usefulness in separating the 
overlapping families. 

A possible improvement in the interloper removal methodology from a 
different perspective was proposed by Radovié et al. (2017). These authors 
developed a variation of the standard HCM approach, introducing an addi- 
tional step in the procedure to reduce the above-mentioned chaining effect. 
The main idea is to prevent chaining through an already identified interloper. 
A family is first identified using the classical HCM approach in this two-step 
procedure. Then, potential interlopers identified based on their reflectance 
characteristics (e.g. spectrum, colour, albedo) are removed from the input cat- 
alogue, and the HCM is rerun. This prevents connecting new members to the 
family by linking them to a suspected interloper. The obtained results show 
that the number of potential interlopers among family members could be sig- 
nificantly reduced in this way. However, a drawback is that this approach could 
miss a part of the family for families of a very complex shape. 

Therefore, information about surface reflectance characteristics helps refine 
the membership of an asteroid family. Still, a sort of caution is needed here. 
Interloper removal methods based on the composition of asteroids associated 
with a family intrinsically assume a homogeneous family. Though available 
data indeed suggest that most families are homogeneous in composition, there 
are some exceptions to this. Such an approach that is apriori giving up the 
possibility of a heterogeneous family could prevent us from discovering hetero- 
geneous families. It contradicts, in some way, claims that families provide us 
information about the parent's body composition. Braking this circular depen- 
dence is not trivial, but a determined fraction of interlopers might be a hint 
if we deal with a heterogeneous family originating from a partly differentiated 
body. Based on the statistical arguments, Migliorini et al. (1995) found that the 
fraction of interlopers in an asteroid family identified using the HCM should 
be up to about 1096. Any fraction significantly larger than that could point 
out to a family heterogeneous in composition, or overlapping of two different 
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families. This should be considered when a level of homogeneity /heterogeneity 
among family members is analysed. 


2.3 Recent advances in asteroid family identification 


Dealing with big data sets has become the most challenging aspect of asteroid 
family identification in recent years. From a technical point of view, the classi- 
cal HCM applied to a large number of asteroids becomes a very time-consuming 
procedure and, at the same time, requires a significant amount of available 
memory. Therefore, the problem is related to the computing resources needed 
to deal with many asteroids. Additionally, a large number of asteroids produce 
a high number-density of asteroids in the space of orbital elements causing 
many families to overlap. This is because the smaller family members are both 
launched at higher ejection velocities and faster transported away from the 
centre of the family due to the non-gravitational effects (see Section 3). 

'There are two main directions proposed so far to cope with the problem. A 
possible solution to deal with a large amount of data for the purpose of classi- 
fication of asteroids into families could be a multi-step methodology proposed 
by Milani et al. (2014, see also Milani et al. (2016)). Alternatively, the prob- 
lem could be treated by employing machine learning-based tools, which can 
be used either to identify new asteroid families (Carruba et al., 2019; Vujicié 
et al., 2020), or to attach new members to known families (Carruba et al., 
2020a). 

Another point of view was brought into the game by Walsh et al. (2013, see 
also Delbó et al. (2017)). Instead of looking into the 3-D space of the proper 
orbital elements, the authors proposed identifying families using the so-called 
V-shape. The method is especially suitable for searching for very old and, 
therefore, highly dispersed families. 

It also becomes increasingly important to find the best way to include 
different datasets available on asteroids. It refers to the development of the 
algorithms able to extract the maximum information from the available data, 
but also treat missing data as the different types of data are available for a 
different number of asteroids. A step forward in using an extended number 
of parameters to identify mutually related objects and treating missing data 
was recently proposed by Holt et al. (2021). The authors employed a cladis- 
tic approach to classifying the Jupiter Trojans. The main advantage of this 
method is that not all characteristics need to be known for an analysis to be 
carried out. This allows using more parameters without truncating the dataset 
due to missing values. 

In the following subsections, we describe the main properties of these new 
strategies and discuss their most important advantages and limitations. 


2.3.1 Multi-step approach 


A multi-step methodology in asteroid families identification is related to par- 
titioning the asteroid data set into subsets typically based on asteroid sizes 
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(absolute magnitudes). The general idea is to apply the HCM to larger objects 
first, and then include the smaller ones in subsequent steps. The most com- 
plete multi-step method was developed by Milani et al. (2014). Their procedure 
consists of the following main steps: 

1. Divide the main belt into several zones to reduce the total number of 
objects which need to be treated simultaneously. 

2. Remove objects having absolute magnitude H fainter than a given 
threshold Hii. The limit typically corresponds to the observational com- 
pleteness limit in a given zone, though in principle, other choices are possible 
as well. 

3. Apply the HCM to each zone to the samples of asteroids with H < Hii. 
In this step, the so-called family cores are identified. 

4. Classify individual asteroids which had not been used in the core clas- 
sification by attaching some of them to the established family cores. In this 
step, the asteroids having a distance from at least one core family member not 
larger than the critical distance are attributed to the core families. 

5. Apply the HCM in each zone to the samples of intermediate background 
asteroids, defined as the set of all the objects not attributed to any family in 
steps 3 and 4. 

6. Attach remaining background asteroids to all families using a single-step 
HCM. In this step, asteroids not associated with any family in steps 3 and 
5 are attributed to all of these families if a distance from at least one family 
member is below the given critical distance. If an asteroid is attributed to more 
than one family, it is considered part of an intersection. 

7. Merging halo families with core families. There are essentially two pos- 
sible cases for families identified at step 5: these families can either be fully 
independent, new families having no relation with the families identified previ- 
ously, or they may be found to overlap with families from steps 3 and 4, forming 
haloes of smaller objects surrounding some family cores. These two cases are 
distinguished based on information about intersection objects identified in step 
6. 

'The above-described procedure is shown to be efficient in dealing with rela- 
tively large datasets. Milani et al. (2014) applied it to about 330,000 numbered 
asteroids, while Milani et al. (2016) extended the analysis to multi-opposition 
objects yielding in total more than 500,000 asteroids. It also solved both prob- 
lems, identification of new families and attribution of new asteroids to existing 
families. 

Still, in its present form, it would be very challenging to successfully apply a 
multi-step method to an order of magnitude more asteroids, which are expected 
to be discovered by ongoing and forthcoming sky-surveys over the next several 
years. A possible solution could be a further partition of the dataset and the 
introduction of new steps in the procedure. That would, however, probably 
make the process more difficult to follow and less reliable. 

'The main problem with this idea is that, in some cases, it splits members of 
an asteroid family into two (or even more) independent groups. Indeed, Milani 
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et al. (2014) suggested a pathway towards clarification of such cases, but it 
could take years of monitoring to resolve such cases. Additionally, merging 
separate families based on the number of intersection asteroids is sometimes 
ambiguous and, at least in part, based on subjective judging. 


2.3.2 Machine learning-based methodology 


In the last decade, machine learning (ML) has emerged as one of the most 
promising complements to traditional data analysis and modelling methods in 
scientific fields. The first steps towards applying ML-based tools for identifi- 
cation of asteroid families and attribution of new members to existing families 
have been performed.? 

Carruba et al. (2019) were the first to attempt to employ ML-based tech- 
niques to identify new asteroid families. The authors investigated the prospect 
of employing ML-based clustering algorithms available in the PYTHON pro- 
gramming language for the identification of the asteroid families. In this work, 
Carruba et al. (2019) adopted the standard distance metric proposed by Zap- 
palà et al. (1990) and focused on a sub-population of main-belt asteroids 
moving onto the high-inclination orbits. This preliminary investigation reveals 
the main advantages and disadvantages of utilizing ML tools for asteroid fam- 
ily identification. On one side, the algorithms are shown to be able to identify 
all significant groupings in the region and to provide results in shorter times 
than the traditional HCM approach. The ML-based HCM turns out to be also 
easier to use than the traditional one. 

On the other hand, regarding the assessment of the membership of indi- 
vidual families, some limitations of the tested ML-based algorithms become 
apparent. The precision of the algorithm, defined as a ratio between the num- 
ber of correctly identified family members and the total number of objects 
associated with a family, varies significantly from family to family. Though pre- 
cision of 10096 has been achieved in some cases, it varies from 37 to 10096, with 
a membership precision below 8096 for more than a half of families (Carruba 
et al., 2019, Table 4). 

Another point to consider is that Carruba et al. (2019) tested ML-based 
algorithms to identify families in less populous, high-inclination regions of the 
main belt. In this region, however, even the traditional HCM is still applicable 
due to the significantly smaller number of objects than the main belt's low- 
inclination part. 

Vujicié et al. (2020) explored the possibility of using an artificial neural 
network (ANN) for the classification of asteroids into families. The obtained 
results showed that the ANN-based algorithm performs somewhat better than 
the clustering algorithms tested by Carruba et al. (2019). In particular, the 
ANN was shown to be efficient in reproducing the traditional HCM results 
for large families with more than 1,000 members. The results for smaller 


?For a review on a more general application of machine learning to asteroid dynamics we refer 
the readers to review by Carruba et al. (2021). 
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families are found to be less reliable, suggesting that ANN may not consis- 
tently outperform the ML-based algorithms, and good practice would be to 
try both approaches whenever possible. Àn additional limitation of an ANN- 
based method is that it acts as a black box and that different runs could lead 
to different results that put the reproducibility of the results into question. 

Carruba et al. (2020a) used ML methods to address another problem 
related to asteroid families, that is, the connection of newly discovered aster- 
oids to existing families. This approach bears some similarities with the 
multi-step method. The basic idea is to identify families using only larger fam- 
ily members and then use the ML-based algorithm to attach new members to 
these families. Indeed, in principle, the second step could be applied to any 
family with known membership. Carruba and co-authors tested nine different 
ML-based classification algorithms. They found that the Extremely Random- 
ized 'Trees method performs the best when attaching the new family members, 
typically achieving high precision between 83 and 98% (Carruba et al., 2020a, 
Table A6). It shows that the ML-based algorithms could successfully attach 
new members to asteroid families. The algorithms perform better when attach- 
ing new members to existing families than when classification is done from 
scratch. We also note that supervised ML-based methods may overlook new 
patterns in the distribution of family members. These algorithms, trained on 
the orbital distributions of family members of known families, will look for 
similar patterns and therefore may not identify unusual structures. 

In Fig. 1 we show how the algorithms for attribution of new family mem- 
bers to an already existing family (ML-Based extra tree and one-step HCM 
algorithms) compare to each other, and how concerning the classical HCM 
approach. Applying the classical HCM at a cut-off of 18 m s^! to an older (and 
therefore smaller) catalogue of proper elements that contains 507,449 objects, 
and to a larger catalogue with 793,310 asteroids, we identified 1452 and 2988 
Veritas family members, respectively. In this case, applying the classical HCM 
to the catalogue that contains about 5696 more asteroids than the initial one, 
and keeping the cut-off distance at 18 m s^, increased the family membership 
for 1536 members, i.e. by a factor of 2. 

On the other hand, applied to the extended catalogue, the Extra 'Tree and 
one-step HCM algorithms respectively attributed 917 and 908 members to the 
family initially identified using the classical HCM and the smaller catalogue. 
From Fig. 1, it is apparent that the distribution of members attributed by the 
Extra Tree (orange squares in the upper panel) and one-step HCM (magenta 
squares in the lower panel) algorithms are almost identical. Therefore, these 
two algorithms show very similar performances. 

The classical HCM algorithm attached a significantly larger number of 
objects when applied to the larger catalogue. However, the distribution of the 
identified members (grey circles in Fig. 1), shows very irregular patterns. As 
discussed in Section 2, this is at the same time related to the main advantage 
and the main disadvantage of the classical HCM. Though many families are 
definitely of irregular shapes, some of these structures might not be the real 
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Fig. 1 Identification of Veritas family members using different methods and different cata- 
logues. In both panels, grey circles represent the family members identified by the classical 
HCM approach applied to the catalogues containing ~793k objects, respectively, and using 
the cut-off of 18 m s~!. Orange (upper panel) and magenta (lower panel) squares represent 
the members attached to the family from the more extensive catalogue, to the members ini- 
tially identified with the classical HCM at 18 m s~! applied to the smaller catalogue with 
~507k records. The upper panel (orange squares) are objects attached by the ExtraTrees 
ML-based algorithm with a hyper-parameter (number of estimators) equal to 41 (Carruba 
et al., 2020a). In the lower panel (magenta squares) are shown members attributed to the 
family with the single 18 m s~! step HCM approach (Milani et al., 2014). 


parts of the families, but rather interlopers associated with the family due to 
the chaining effect. The most severe problem that may arise due to chaining 
is merging two or more separate families. 

Given the above example, both ML-base extra tree and one-step HCM 
algorithms perform pretty well in attaching new members to an existing family. 
Therefore, we recommend using one of these algorithms for fast checking of 
new family members. The algorithms could be especially useful in the areas 
of a very high number-density of asteroids. On the other hand, the classical 
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HCM approach could still overperform the two attribution algorithms. It still 
seems like the best choice in cases where the number-density of asteroids and 
distances between nearby families allow its application. 


2.3.3 Primordial families and V-shape based detection of 
families 


The classical HCM approach typically fails to identify very old families, formed 
more than about 3 Gyr ago. Such ancient groups are so dispersed in the space 
of proper orbital elements that the HCM cannot distinguish them from the 
background population of asteroids. For this reason, many current efforts have 
focused on detecting ancient families using different techniques. 

Recent methods make use of the V-shape (see Section 3.2) to identify fam- 
ilies that are too dispersed to be discovered using the classical HCM. The idea 
is to test for possible V-shape number-densities amongst the background of 
unassociated asteroids. V-shape finding algorithms can be divided into two 
categories. The first is the so-called border method (aka dC-method or dK- 
method) that utilises a ratio between the number of objects just outside and 
inside the V-shape. It is helpful when families have a distinct edge and stand 
out with respect to the background asteroids (see Walsh et al., 2013; Delbó 
et al., 2017). The second category refers to the density methods (aka aw- 
method). The goal is to detect a peak of asteroid density in a V-shape plane. 
It looks like a better technique for families embedded in a dense background 
of asteroids (Bolin et al., 2017). 

These V-shape based family identification techniques are not trivial to 
implement. Here we review the main steps, while for full implementation 
details, we refer readers to Bolin et al. (2017) for the implementation of the 
dC, and Delbó et al. (2017) for the a, version of the method. 

The dC-method or dK-method is part of the border method category. It is 
based on testing a family centre and age or slope (how large the V-shape is). 
This technique tests a set of family centres and family slopes, evaluating a ratio 
between the bodies just inside and outside the defined V-shape borders. They 
are called K or C methods because the slope K relates to the parameter C 
through K = \/py/(DoC), where p, is the geometric visible albedo, Do = 1329 
km (Walsh et al., 2013; Delbó et al., 2017), and C is the parameter defining 
the re-scaling asteroid's distance from each other (Vokrouhlicky et al., 2006c). 
The a,-method considers a slightly different approach. It scans for the V- 
shape vertically in the (a5, 1/D) or (ay, H) space and calculates the ratio of 
asteroids above and below the tested V-shape. The offset is established as a 
delta in 1/D and named as a, (Delbó et al., 2017). Deienno et al. (2021) 
summarizes the two techniques saying that the dC/dK-method searches for 
(K + AK)|a, — dp e] and the a,,-method searches for (K los, — Ap e] + aw). 

Both categories have been tested on synthetic families (Walsh et al., 2013; 
Bolin et al., 2017; Deienno et al., 2021) and they have also been used to try 
to detect primordial families. Walsh et al. (2013) were the first to use the V- 
shape based family identification finding the new Polana family, about 2 Gyr 
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old group, and the Eulalia family, about 0.8 — 1 Gyr old group, both within 
the Nysa-Polana complex. Delbó et al. (2017) proposed a 4-billion years old 
family extending across the entire inner part of the main belt, whose members 
include most of the dark asteroids previously unlinked to families. Delbó et al. 
(2019) extended the research of ancient families, suggesting other two ancient 
families, Athor and Zita. 

A similar technique was used by Dykhuis and Greenberg (2015), extend- 
ing the work done by Walsh et al. (2013), to study the Nysa-Polana complex 
of families. In addition to seeking clusters in the dynamical and reflectance 
parameter space, the authors seek signatures of dynamical evolution that dis- 
tinguish multiple families within a single group. Dykhuis and Greenberg (2015) 
have identified five distinct clusters of asteroids within the Nysa-Polana orbital 
region of the main belt. Two of these clusters are associated with collisions 
on (135) Hertha, two others with (495) Eulalia, and one is linked with (142) 
Polana. 

These different techniques allow identifying some very old families in the 
inner main asteroid belt. However, the V-shape family identification method is 
also very sensitive to many parameters involved in the search, and additional 
verifications are needed to fully quantify the robustness of the method. 

Deienno et al. (2021) have done a comprehensive study to characterise the 
efficiency of the V-shape method in detecting asteroid families. They worked 
using a synthetic population of asteroid families and background asteroids, 
letting those families evolve over billions of years, setting different levels of 
ratio of the synthetic family and background asteroids to derive a detection 
efficiency map for the V-shape method. They also tested their methods using 
the families found by Delbó et al. (2017, 2019). The paper results show that 
families older than 3 Gyr are difficult to detect, while families 0.5-2.5 Gyr 
old are more easily detected, with an efficiency of more than 8096. In light of 
these results, Delbó et al. (2017, 2019) seem to be lucky to detect the ancient 
families. As stated by Deienno et al. (2021), this is still a work in progress, 
even though promising. Other steps will include adding statistical analysis 
and other debiasing techniques to fully characterise the V-shape methods' 
efficiency. Especially, the methods need to be tested for their reliability, that 
is, how often they could potentially generate false-positive results. 


2.3.4 Cladistic approach 


A step forward in using an extended number of parameters to identify mutually 
related objects and treating missing data was recently proposed by Holt et al. 
(2021). The authors employed a cladistic approach to classifying the Jupiter 
Trojans. 

The main advantage of this method is that not all characteristics need to 
be known for a cladistical analysis to be carried out. This allows using more 
parameters without truncating the dataset due to missing values. 

The method, however, has some setbacks as well. First, the algorithm is 
time-consuming and scales badly with the number of objects. Therefore, its 


Springer Nature 2021 IATEX template 


Article Title 13 


extension to the main belt region would be very challenging. Second, though 
the ability to treat the missing data is one of the method's main advantages, it 
should be noted that such asymmetric datasets could lead to spurious or unre- 
liable results in some cases. For instance, if only basic information is available, 
such as apparent magnitudes and orbit solutions, the amount of data is insuf- 
ficient to place the significant weight on the results. Nevertheless, the cladistic 
approach employed by Holt et al. (2021) represents an important step forward 
in analysing mutual relationships among the objects in an asteroid population. 
In order to take full advantage of the approach, caution is needed when inter- 
preting the data. A possible way to increase the robustness of the results could 
be achieved by combining cladistic techniques with other available methods. 


2.4 Asteroid Families Portal 


For the reasons explained above, identifying new families and the associa- 
tion of members to existing families are demanding procedures that require 
careful analysis, and the results of such studies are typically published as 
journal papers. However, a fast-growing number of newly discovered aster- 
oids facilitates also need to constantly update the list of know families, or at 
least to associate new members to existing families. In recent years, web por- 
tals and databases have been recognized as a possible way to cope with the 
large size and complexity of data accumulated over the years. Such portals 
are often devoted to specific purposes and allow scientists to quickly access 
and review different data. With the possibility of integrating the data and 
the computation, such portals enable interactive analytics, typically based on 
well-established tools and algorithms. They, on one side, allow better exploita- 
tion of the available data, while on the other side, make using different tools 
simpler and available to a broader community. 

In 2017, an internet portal on asteroid families was launched, namely 
the Asteroid Families Portal (AFP). The portal is freely available at http:// 
asteroids.matf.bg.ac.rs/fam/, and provides data and tools to study asteroid 
families, including some tools described in Section 2. The primary goal of the 
AFP is to implement various algorithms used to study asteroid families and, 
therefore, to allow users to employ these tools online. Additionally, the portal 
collects different data about families and closely related subjects. 

Currently available services and data at the AFP covers the following five 
categories: 


e Application of the Hierarchical Clustering Method to identify members of an 
asteroid family: as explained in Section 2, the HCM is the most widely used 
method for the identification of asteroid families and their memberships. 
The AFP allows the algorithm described in Zappalà et al. (1990) to be 
run online using either default parameters (for less experienced users) or 
advanced settings (for more experienced users). 

* Application of an automatic procedure to obtain a list of family interlopers: 
this option is based on the algorithm described in (Radovié et al., 2017). It 
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provides the list of potential interlopers and the information used to suggest 
such status of an asteroid. 

e Application of the Backward Integration Method to estimate the age of 
an asteroid family: for a given list of asteroids, the algorithm described by 
Nesvorny et al. (2003) is run over the time series of mean orbital elements. It 
returns, as a result, the time evolution of average differences in secular angles 
(longitude of node Q and longitude of perihelion w) for selected objects. 
'The analysis covers an interval of 10 Myr in the past. 

e Proper orbital elements of numbered and multi-opposition asteroids, and 
active asteroids: the synthetic proper elements available at the AFP are com- 
puted numerically following, in principle, the method described in Knezevié 
and Milani (2003). The important difference with respect to the procedure 
used by Knezevié and Milani, is that orbits of all asteroids are numerically 
propagated for 10 Myr, and using the same dynamical model that includes 
seven major planets (excluding Mercury). Fig. 2 shows the snapshot of the 
page at the AFP devoted to the proper elements. 

e List of the peer-reviewed papers on asteroid families: the currently available 
data goes back to 2007, intending to extend it more in the past. 


of January 2021) 


the specified values: (19.2.12 7,2111.7,1126,2013.6,111,147) 
fed values: () 


sed to read the specified values: (22.9.2127 2111.7.112.6,2113.6,/11.1,7) 
read the specified values: (a22,13.7,2112.7,112.6,2013.0) 


Fig. 2 Snapshot of the Proper Elements page at the Asteroid Families Portal. The page is 
available at: http: //asteroids.matf.bg.ac.rs/fam/properelements.php 


2.5 Future prospects and challenges for asteroid families 
identification 


Identifying families and classifying individual asteroids into families is already 
a very challenging task. In a way, paradoxically, these challenges are mainly 
due to a large amount of newly available data, especially newly discovered 
asteroids. The problems will soon become even additionally severe once more 
extensive data sets produced by the large sky-surveys become available, espe- 
cially when the Vera Rubin telescope will start observing (Ivezić et al., 2019; 
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Vera C. Rubin Observatory LSST Solar System Science Collaboration et al., 
2021), discovering millions of new small solar system objects. 

On a purely computational side, the classical HCM method is not suitable 
for dealing with many objects. It is primarily not due to long computation 
time but rather because of extreme memory requirements. These limitations 
are challenging, but there are opportunities to overcome them, especially the 
execution times, by newly available clustering techniques. Even already devel- 
oped methods presented in Section 2.3 could overcome these issues to some 
degree. 

A more difficult problem is a large number density of asteroids, which 
reduces the average distance between two objects below the uncertainty in 
the distance computation. A large number of objects will also further increase 
the chaining effect, which will cause families to overlap in the space of orbital 
elements. These problems will be tricky to solve. 

Despite being typically of much better relative accuracy than the phys- 
ical parameters on asteroids, uncertainties in the proper elements further 
complicate the situation. In Fig. 3, we show the distribution of the position 
uncertainties caused by uncertainties in proper orbital elements. The uncer- 
tainty in the position of an individual asteroid is computed using the standard 
HCM metrics given by Eq. (1). These uncertainties in the position of individual 
objects translate into uncertainties in the distances between the asteroids. If 
significant, such uncertainties further complicate classifying asteroids into fam- 
ilies. Until recently, these uncertainties were safely below the typical distances 
among the asteroids in the 3D space of proper elements. However, a recent 
increase in the number of known asteroids significantly decreased the average 
distance between the two neighbouring objects, making position uncertainty 
an important issue. It needs to be taken into account in future works related 
to identifying asteroid families and attributing new members to known groups. 

Another aspect of the family identification problem is the generally 
unknown distribution of the background (non-family) objects. A standard 
assumption behind applying clustering algorithms to identify asteroid fami- 
lies is the uniform distribution of the background (non-family) objects. It is 
well known, however, that this is not the case. There are at least two rea- 
sons for that. First, dynamical evolution removes faster objects from some 
regions of the main asteroid belt than from others. Second, the procedures 
used to compute the synthetic proper elements introduce some non-uniformity 
by themselves. Good examples are concentrations of asteroids near the centres 
of the resonances, which is a consequence of an averaging procedure used to 
compute the proper elements. A potential way to reduce the effect of this issue 
is to use a synthetic background that better represents a random population 
in the proper elements space. 

Such a population could be constructed starting from a uniform distri- 
bution in osculating elements space and computing the corresponding proper 
elements for a number of test particles. The outcome of such an approach is 
shown in Fig. 4. 
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Fig. 3 Position uncertainties distribution in the space of the proper orbital elements, 
expressed in meters per second. The data for about 507,000 numbered asteroids are shown. 
The errors are computed according to Eq. (1), using the formal uncertainties in proper 
semi-major axis, eccentricity and sine of inclination. 
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Fig. 4 Proper elements distribution for about 200,000 test particles uniformly distributed 
across the main belt in the osculating orbital elements. 


The most difficult to solve is probably the problem of the overlapping 
families. To deal with this issue, we need to look in another direction. For 
instance, identifying families in multi-dimensional space (N73) could be an 
option, but appropriate variables and metrics need to be defined. For instance, 
combining the standard metrics with the V-shape criterion could help a bit. 
However, this would apply only to families older than about 50 Myr. Also, 
for one-side families, it may not work well. On the other hand, in the case 
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of overlapping families, the main problem is not necessary to recognize that 
there are more than a single family, but the issue is how to attribute individual 
asteroids to these groups correctly. These remain open challenges for future 
works. 

Certainly, identifying asteroid families and attribution of new members to 
existing families will be an important task in the near future. We discussed 
here the main problems and proposed some of the directions towards their 
solution. However, more comprehensive work is needed to understand better 
which techniques and algorithms will or will not work. 


3 Long-term dynamical evolution 


'To successfully reconstruct the long-term dynamical evolution of an asteroid 
family, it is essential to understand all the relevant mechanisms that change the 
family and introduce these effects into a model. The relevant dynamical mech- 
anisms could be broadly divided into the gravitational and non-gravitational 
effects. 'To account for the gravitation effects, generally, the best strategy would 
be to include the gravitational influences of all the solar system objects in the 
model. Due to computational reasons, this is, however, not feasible. Therefore, 
a balance between the approach needs to be followed. Luckily, for the main-belt 
asteroids, gravitational perturbations of most of the solar system objects are 
negligible. The relevant perturbations are typically caused by eight major plan- 
ets and a few of the most massive asteroids. In some cases, even not all these 
objects need to be taken into account. Regarding the non-gravitational effects, 
for asteroids, relevant are the Yarkovsky and Yarkovsky-O'Keefe-Radzievskii- 
Paddack (YORP) effects. Except due to the dynamical evolution, families 
evolve for other reasons, and collisional evolution is especially important. 
'The dynamical and collisional evolutions are coupled, and simulations of the 
dynamical evolution often need to account also for the collisional evolution. 

'The simulations of the long-term evolution of asteroid families have seen 
many remarkable achievements over the years. There is a long list of literature 
covering these topics. It is beyond our goal to mention all these works. Instead, 
we primarily focus on the results published since 2015 and briefly describe 
some of the key results from the more distant past. A more in-depth review 
of the results achieved in the more distant past can be found elsewhere (see 
for instance Cellino et al., 2009; Nesvorny et al., 2015, and references therein). 
For a review of the dynamics of families interacting with secular resonances, 
we refer the readers to the recent review by Carruba et al. (2018a, see also 
Knezevié (2022)). 

Probably the most important step forward in modeling the long-term orbit 
evolution of asteroid families, made in the last several years, was an intro- 
duction of the stochastic YORP concept. The concept, based on the work by 
Statler (2009), was introduced by Bottke et al. (2015) in their search for a pos- 
sible origin of asteroid Bennu. The authors used a Monte Carlo approach to 
study the dynamical spreading in terms of the orbital semi-major axis of five 
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inner main-belt families. Bottke et al. (2015), however, assumed that due to 
shape changes to asteroids, produced by processes such as crater formation or 
changes to asteroid rotational angular momentum by YORP, may cause the 
asteroids’ spin rates to undergo a random walk. The stochastic YORP mech- 
anism slows down how often asteroids reach YORP end-states. The effect, in 
turn, allows the faster drift of the asteroids' semi-major axis than in the clas- 
sical Yarkovsky/YORP scenario. This new model successfully reproduced the 
semi-major axis distribution of the observed families, and opened possibilities 
to explain some other mismatching related to the modeling of the YORP effect. 

Among the recent highlights was also the discovery of the dynamical impor- 
tance of secular resonances with massive bodies other than the major planets. 
Novakovié et al. (2015) found that the spreading in orbital inclination seen in 
the Hoffmeister family is an outcome of the vic = s — sc nodal linear secu- 
lar resonance with Ceres. The authors showed that passing through the vic 
resonance may cause significant changes in the orbital inclination of an aster- 
oid. It was the first direct proof that a secular resonance between Ceres and 
other asteroids, previously completely overlooked, can cause significant orbital 
evolution. 

'To illustrate the principle of the dynamic evolution simulation and better 
explain the role of secular resonances with Ceres, we discuss the latter in more 
detail. For this purpose, some possible scenarios of the dynamic evolution of 
the Hoffmeister family members are shown in Fig. 5. The figure shows the orbit 
evolution of test particles over 200 Myr. Three representative test particles are 
selected with negative Yarkovsky drift and, therefore, drifting toward the Sun. 
They are selected to show the interaction with the vic secular resonance. In 
the top panels, we see an example of the particle drifting only in terms of the 
proper semi-major axis for about 65 Myr, when it reaches the vic resonance. 
'The interaction with the resonance increased the orbital inclination of the test 
particle. Once the particle exits vic, the Yarkovsky effect drives it to zı in 
about 60 Myr. After that, zı offers a path followed by the asteroid until the 
end of the simulation. 

'The middle panels show the evolution of an asteroid with the fastest neg- 
ative Yarkovsky drift used in the simulations. In this case, the particle went 
across the resonance in a relatively short time, about 10 Myr, and interaction 
with the vic only slightly increased the orbital inclination. Once outside the 
resonance, the particle continues drifting in the semi-major axis. The orbital 
eccentricity and inclination are not affected, suggesting that the particle is not 
interacting with the zı resonance. 

'The bottom panels illustrate the evolution of an asteroid whose proper incli- 
nation is significantly reduced by the uc resonance. It reduces approximately 
0.021 in the sine of proper inclination, meaning a loss of 1.2 degrees in proper 
inclination in less than 100 Myr. Such a scenario is favoured by the orientation 
of the 1c, which creates a path in agreement with the inward direction of the 
Yarkovsky drift. While the inclination is dropping down, eccentricity remains 
approximately constant. This behaviour is in agreement with the fact that the 
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Fig. 5 Three representative scenarios (top, middle and bottom panels) of the orbit evolu- 
tion of an asteroid from the Hoffmeister family. The left and right panels show the evolution 
in the semi-major axis vs inclination (ap, sin(ip)) and the semi-major axis vs eccentricity ap, 
ep) planes, respectively. The examples highlight the interaction of the Yarkovsky effect and 
the vic secular resonance with Ceres. The dashed dark red line represents the approximate 
location of vic resonance. The full and dashed blue lines represent the centre and one border 
of z1 secular resonance, respectively. Thin vertical dashed lines indicate the mean motion 
resonances. Grey dots represent the distribution of the family members after 200 Myr evolu- 
tion. Orange squares show the successive locations every 10 Myr of a selected test particle. 
When the vic catches the particle, the orange squares change into dark red triangles, while 
inside the z1, the particle is represented by blue circles. 


particle remains far enough from z1, so that this resonance cannot influence 
the motion of the particle. 

The work by Novakovié et al. (2015) opened several new possibilities to 
investigate small bodies! dynamics and initiated related studies. An interesting 
aspect in this respect is that by modifying the orbital inclination of family 
members, the nodal secular resonances indirectly modify the initial ejection 
velocity field, particularly an out of the orbital plane velocity component. As a 
result, the distribution of this velocity component for members of an asteroid 
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family becomes more and more leptokurtic? over time (Carruba et al., 2016c). 
This fact could be, in turn, used to estimate the initial (post-impact) velocity 
distribution. For details on the studies along these lines, see Section 3.1. 

The Astrid family members have an unusual orbital inclination distribu- 
tion, the feature caused by an interaction of the family with the vc nodal 
secular resonance with Ceres (Novaković et al., 2016). Consequently, the 
observed distribution of the vy component of terminal ejection velocities is 
quite leptokurtic. Carruba (2016) used the orbital configuration of the Astrid 
family to set constraints on the parameters describing the strength of the 
Yarkovsky force, such as the bulk and surface density and the thermal conduc- 
tivity of surface material. By varying the initial conditions, and by demanding 
that the current value of the kurtosis of the distribution in vy be reached over 
the estimated lifetime of the family, Carruba (2016) found the thermal con- 
ductivity of the Astrid family members should be 0.001 W m7! K-t, and that 
the surface and bulk density should be higher than 1000 kg m”*. Furthermore, 
using the Monte Carlo approach to simulate the Yarkovsky/YORP evolution 
of the family, the author estimated its age to be 140 + 30 Myr. Similar argu- 
ments were used by (Carruba et al., 2017a) to reassess the Hofmeister family. 
By combining the results from a Yarkovsky/Y ORP Monte Carlo method with 
the constraints provided by the inclination distribution and the time behaviour 
of the kurtosis of the vy component of the ejection velocity field, the authors 
found that the Hoffmeister family should be 220*% Myr old. 

Carruba and Barletta (2019) studied the dynamics of the Nemesis asteroid 
family, also considering the influence of Ceres. An especially unusual character- 
istic of the family is the asymmetric distribution of its members in terms of the 
orbital semi-major axis. More than 8596 of the members are located at smaller 
semi-major axes with respect to the centre of the family. The authors, however, 
found that the asymmetric shape of this family is probably an outcome of the 
original ejection velocity field, rather than the post-impact dynamical evolu- 
tion. Considering the population of family members currently present in the 
Vic resonance, Carruba and Barletta (2019) set an upper age of the Nemesis 
family to be 200 Myr. 

Regarding the interaction with vic secular resonance, an interesting case 
for further study is the Seinajoki asteroid family. The family is also crossed 
by this resonance (Tsirvoulis and Novakovié, 2016), and its overall shape in 
the (a,,sinip) plane could be explained by the effects of the uc resonance 
(Novakovié et al., 2016). However, although the perturbation caused by this 
resonance explains the overall shape of the Seinajoki, several smaller features 
remain unexplained. Spoto et al. (2015) found that the V-shape boundaries of 
the family are consistent with two different ages, suggesting that the family 
might be an outcome of two collision events. 


3A leptokurtic distribution has positive kurtosis larger than that of a corresponding Gaussian 
distribution. In general, leptokurtic distributions have heavier tails or a higher probability of 
extreme outliers values when compared to Gaussian distributions. 
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The proper orbital elements may be especially suitable for studies of 
the long-term dynamical evolution of asteroid families. These are the quasi- 
integrals of the full N-body equations of motions and, as such, nearly constant 
in time, thus preserving the dynamical signatures from the remote past (see 
e.g. Knezevié, 2017, and references therein). For these reasons, the proper ele- 
ments (or their approximation as long-term averages) are often used to follow 
the dynamical evolution of asteroid families. 

Nesvorny et al. (2002a) were the first to study the long-term evolution of 
an asteroid family in the space of proper elements. The authors used averaged 
elements over a long time as a proxy for the proper orbital elements. The long- 
term evolution of the Flora family was simulated by combining shorter N-body 
simulations with a kind of Monte Carlo simulation. In this paper, Nesvorny 
et al. focused on the dynamical dispersion of the proper eccentricity and incli- 
nation, which occurs due to the chaotic diffusion in narrow mean motion 
resonances, and considered the Yarkovsky non-gravitational force and the 
gravitational impulses received at close approaches with large asteroids. More 
recently, Novakovié et al. (2010b) employed a similar Monte Carlo approach 
to study the evolution of the Lixiaohua family in the space of proper elements. 
The authors used a model which combines the chaotic diffusion, Yarkovsky 
thermal effect, and random walk due to close encounters with nearby massive 
asteroids. These simulations showed that all the effects should be taken into 
account to reproduce the observed distribution of family members accurately. 

Still, although the two above analyses were performed mimicking the evo- 
lution in the space of proper elements, these were based on Monte Carlo 
techniques. The evolution of asteroid families in the space of proper elements 
but based on the full N-body simulations were performed by Novakovié et al. 
(2015) to study the Hoffmeister family, and more recently by Novaković et al. 
(2017) to study the Tamara family. 

Recently, Broz and Morbidelli (2019) developed an improved method to 
study the 3-dimensional shapes of asteroid families in the space of the proper 
elements. The method relies upon comparing the observed and synthetic fam- 
ilies from N-body simulations as directly as possible. A primal novelty of the 
approach is in simultaneous modelling of the evolution of both a family and 
background population of asteroids. Broz and Morbidelli (2019) applied their 
model to the Eos family and were able to explain most of the family's features 
fully. 

The long-term evolution of asteroid families due to non-gravitational effects 
may produce specific patterns in family members’ distribution. The well-known 
phenomenon is the formation of family ears (Vokrouhlicky et al., 2006a; Lowry 
et al., 2020). The ears form because the small family members concentrate 
toward the extreme values of the semi-major axis. At the same time, there is a 
lack of asteroids in the family center, giving the family the appearance of ears. 
'This feature is a consequence of the YORP effect, which produces a perpen- 
dicular spin axis orientation relative to the orbital plane and maximizes the 
Yarkovsky drift. T'his expected central depletion is, however, not always clearly 
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noticeable. Paolicchi and Knezevié (2016) investigated this feature in detail 
and developed a concept of the YORP-eye and a central depletion parameter. 
The analysis allowed the authors to find YORP footprints in several families. 
Furthermore, Paolicchi and Knezevié (2016) also found that the location of 
the YORP-eye within the family evolves with a family. Therefore, the YORP- 
eye could be used as an indication of the age of the family (see Section 3.2 for 
further details on the age estimation methodologies). 

An excellent example of what we can currently achieve in modelling the 
dynamical evolution of an asteroid family is a recent study of the Clarissa 
family by Lowry et al. (2020). In this work, the authors primarily aimed to 
apply the Yarkovsky effect chronology to estimate the age of the family. In 
other words, the goal was to define the appropriate initial conditions and follow 
the dynamical evolution in time until the observed orbital structure of the 
Clarissa family is matched. The unique challenge, in this case, was that it is 
highly asymmetrical in the proper semi-major axis. Therefore, to successfully 
reproduce the shape of the family, Lowry et al. ran multiple simulations by 
varying initial ejection velocity fields and the fraction of initially retrograde 
spins. 

Additionally, the authors tested different Yarkovsky/Y ORP models. In par- 
ticular, different YORP models were used by changing how often the effect 
tends to accelerate or slow down the rotation. Interestingly, the best fits were 
obtained with the model assuming a 4:1 preference for spin-up by YORP, and 
assuming that about 8096 of small family members initially had retrograde 
rotation. The corresponding age of the Clarissa family was found to be 56 +6 
Myr. The work by Lowry et al. (2020) demonstrated how a relatively complex 
situation could be solved with careful modelling and extensive testing. 


Limitations of the current models.- One of the most critical challenges that 
we face in modelling the long-term dynamical evolution of families is still the 
proper modelling of the YORP effect. It changes the spin axis orientation and 
spin rate of asteroids and therefore affects the Yarkovsky effect because the 
magnitude of the latter depends on the rotation state of asteroids. However, the 
YORP effect depends on asteroid shapes and small scales topography, such as 
the presence of boulders and craters at the surface. This already complex situa- 
tion is further complicated by the fact that non-destructive collisions could also 
change spin states. It implies that to model the Yarkovsky effect adequately, 
it is necessary to model the combined effect of the YORP and impacts. On 
one side, we are dealing with impacts that are stochastic events and could be 
treated only in a statistical sense. On the other hand, the YORP effect could 
be decomposed into a stochastic and deterministic component. The modelling 
of the deterministic component is improving over the years, and new models 
are developed (see Vokrouhlicky et al., 2015; Golubov and Scheeres, 2019; Gol- 
ubov et al., 2021, and references therein). On the other hand, the contribution 
of boulders, craters, and shape evolution to the YORP effect is generally not 
deterministic. The distributions of craters and boulders, and their evolution 
at the surface of an asteroid are generally unknown. 
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If we only consider the deterministic component of the YORP effect, the 
spin axis evolution due to the YORP typically reduces the magnitude of the 
Yarkovsky effect. It is because, for about a half of asteroids, YORP should 
bring obliquities to 90 degrees (Golubov et al., 2021); see, however, also older 
results of Capek and Vokrouhlicky (2004) which predict the asymptotic YORP 
obliquities near extreme values of 0 and 180 degrees. T'he pole orientation near 
the orbital plane switches off the diurnal component of the Yarkovsky effect, 
which scales with cosy, with y being an asteroid obliquity. For another half 
of objects, the y is pushed towards O or 180 degrees what should, in princi- 
ple, increase the Yarkovsky effect. However, the obliquities are not expected to 
stay at these asymptotic values, but instead, undergo continuous cycles over 
relatively short timescales. As a result, according to this theory, for most aster- 
oids, YORP should reduce the Yarkovsky effect over a time scales significantly 
longer than the YORP cycle timescale. 

Such a scenario, however, does not match the observations. Though not 
related to the families, an excellent example of this is the influx of near-Earth 
asteroids, which is far too low to match the observations when the YORP 
effect is taken into account (Granvik et al., 2017). Therefore, the stochastic 
contribution of the YORP effect and spin evolution due to the impacts need 
to be added to the model. This, however, is not so simple, and authors are 
using different models that yield different results and sometimes even opposite 
conclusions. 

For instance, despite its success in modelling some asteroid families and 
providing a reasonably good match to the observations, the stochastic YORP 
model by Bottke et al. (2015) involves some free parameters which are spe- 
cially designed to justify the model to the observations. These are cyorp 
which defines the strength of the YORP effect, and Creorient that is an empiric 
multiplication factor that controls the characteristic timescale of reorienta- 
tion events. Therefore, the approach is still missing a complete physical 
interpretation. 

Recently, Marzari et al. (2020) presented in some sense the complete model 
of the evolution of asteroid families under Yarkovsky/YORP effects and colli- 
sions. Apart from being the most complete so far, the model has the advantage 
of not involving nonphysical free parameters. However, despite being promis- 
ing at first, the model is yet to be verified against the observational constraints, 
as the authors tested their model only using the synthetic families. However, 
some of the presented results suggest that the observed spreading of the syn- 
thetic families could be far too small to match the spreading of real families, 
except assuming much older ages of these families. For instance, the maxi- 
mum spreading of test particles assumed to be about 5 km in diameter, after 
500 Myr of a Koronis-like synthetic family evolution is only about 0.02 au 
at the inner and even less at the outer side of the family (Fig.10; Marzari 
et al., 2020). This translates into a maximum drift of 2 x 107^ au Myr”!, for 


“See Fenucci and Novakovié (2022) for a recent implementation of these effects into OrbFit and 
Mercury numerical integration packages. 
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a 1 km size body. It is significantly below the theoretically predicted maxi- 
mum Yarkovsky induced drifts or those derived, for instance, by scaling from 
the value obtained for asteroid Bennu, which suggests a maximum drift faster 
than 3 x 1074 au Myr”!. Indeed, Marzari et al. (2020) used the combined 
Yarkovsky/YORP model, which is expected to produce a somewhat smaller 
semi-major axis drift than the purely Yarkovsky effect. However, the level of 
the slowdown seems not to be supported by observational evidence. 

Also, Marzari et al. did not observe any formation of the YORP-eye. They 
suggested that this feature could be related to collisional physics rather than 
to the subsequent evolution driven by the Yarkovsky effect. However, this fact 
could also point out the model's limitations. 

Despite all the above-described efforts and remarkable results in recon- 
structing the past evolution of asteroid families and their members, there are 
still many open questions and models to be improved. At the same time, 
it is essential to realize that asteroids have been subject to a large number 
of diverse stochastic events. Therefore, subsequent collisional and dynamical 
processes have masked information about the fragments produced by past col- 
lisions. It means the initial conditions and the post-impact evolution for many 
asteroid families (and especially the old ones) may never be precisely known. 
Given these limitations, modellers do the best they can with available informa- 
tion. Choosing parameters and formalism within the bounds of what is known 
and testing the results against the available constraints is a recipe to follow. 
Though, the interpretation of even good matches must always be met with 
some scepticism and wariness. Besides, a careful modeller needs also to run 
simulations over numerous trials to characterize how random events may have 
influenced outcomes. 


3.1 Constraints on the initial ejection velocity fields 


Reconstruction of an initial ejection velocity field (EVF) of asteroid families is 
essential to better understand the collisional events’ physics. Such information 
could be cross-validated with the outputs from impact hydrocodes simula- 
tions (Michel et al., 2015). In turn, this may help, for instance, improve the 
asteroid internal structure models. However, evaluating an initial EVF of an 
asteroid family is a complex issue. Having said that, let us note that the EVF 
represents an initial state in the simulations of the long-term dynamical evo- 
lution. A reasonable assumption on the size of EVF could be made based on 
the results from the numerical simulations of asteroid disruptions. Still, more 
accurate data for a specific family needs to be extracted from what we know 
about this family. However, the EVFs are often determined from the long-term 
evolution simulations, similar to the family ages. Therefore their estimation 
is coupled with both the evolution process and age determination. It implies 
that, in principle, family age and EVF could be simultaneously determined by 
simulating the evolution of the family until it matches the available observa- 
tional constraints. T'his is, however, tricky, and much better results could be 
obtained if age is somehow independently constrained. 
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'Therefore, while modelling a post-impact evolution of an asteroid family 
allows, in essence, reconstruction of the family's initial ejection velocity field, 
its efficiency depends from case to case. This approach works the best for young 
asteroid families, though in some cases, reliable results could also be obtained 
for older families. 

A pioneering attempt in defining the size-ejection velocity relationship for 
asteroid families! members was made by Cellino et al. (1999). The authors 
found an inverse dependence of ejection velocities on the size, characterised by 
an exponent between —2/3 and —1. While successful in explaining the avail- 
able data on asteroid families at that time, its important limitation was that 
it did not account for the family evolution. The result could be at least in 
part compromised by the fact that the Yarkovsky effect is causing a similar 
size-dependent evolution. Therefore, the identified size-ejection velocity rela- 
tionships might be consequences of the post-impact evolution due to this effect 
rather than an outcome of the initial collision event. Even more sensitive to 
this analysis are the obtained ejection velocities, which may be overestimated 
by Cellino et al. (1999). However, we note that this analysis was done using 
objects larger than about 10 km in size, a population known at that time. 
These are large enough objects that, in some cases, their Yarkovsky-induced 
post-impact evolution should be relatively limited. T'herefore, the result could 
be still valid to some degree. 

Among the subsequent work, let us mention a series of papers by Vokrouh- 
licky et al. (2006b,c,a). These authors reconstructed the post-impact evolution 
of several asteroid families and found that the initial EVFs were significantly 
smaller compared to the present-day situation. In particular, a typical ejec- 
tion velocity of 5 km fragments in the Erigone, Massalia, Merxia, Astrid, and 
Agnia families was found a few tens of meters per second. This value means 
the initial semi-major axis dispersal in analysed families was 25 — 5096 of the 
observed value. The post-impact evolution naturally explained such a situa- 
tion, and newly derived initial velocities match well the results of hydrocode 
simulations and laboratory experiments (Michel et al., 2015). 

For results derived from young asteroid families, we refer interested readers 
to Sections 4.1 and 4.2. However, let us recall here the results obtained by 
Nesvorny et al. (2006a) for the Karin asteroid family as those changed the 
paradigm of how large the initial velocity fields of asteroid families. T'he Karin 
was the first young asteroid family discovered. Being just 5.8 Myr old, the 
evolutionary mechanisms did not have enough time to modify Karin family 
fragments’ velocity field significantly. Therefore, the members of the Karin 
family provide direct insight into its original EVF. 

Nesvorny et al. (2006a) performed a set of hydrocode simulations of the 
Karin family formation, and found good fits to the size-frequency distribution 
of the observed fragments and the ejection speeds inferred from their orbits. 
The results suggested that the ejection speeds of smaller fragments produced 
by the collision were faster than those of the larger fragments, with the mean 
values for > 3-km-diameter fragments being about 10 m s”!. 
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Despite this success, it was clear that precise determination of the initial 
EVFs is a complex issue and that, especially for older families, novel method- 
ologies are needed to deal with these challenges. The problem was recently 
reassessed by Carruba and Nesvorny (2016). The authors investigated the pos- 
sibility to use only the vertical component vw of the EVF, which is related 
to the orbital inclination. The motivation behind this approach is that the 
proper orbital inclination is generally less affected by the dynamical evolu- 
tion processes than the orbital semi-major axis and eccentricity. Carruba and 
Nesvorny (2016) found that the distribution of the vw should be more peaked 
than a Gaussian distribution (i.e. be leptokurtic) even if the initial distribu- 
tion were Gaussian. They surveyed known asteroid families to distinguish such 
cases and identified eight families characterized by the leptokurtic vw distri- 
bution. All these families are located in dynamically quiet regions of the main 
belt, and their initial orbital configuration should almost not be modified by 
subsequent orbital evolution. Therefore, these families are good candidates for 
the reconstruction of the initial EVF. 

One of these eight identified families is the Koronis family, and in the follow- 
up work, Carruba et al. (2016a) investigated the initial velocities in this group. 
'The main results are that the spread in the original ejection speeds is consis- 
tent with an inverse size-velocity interdependence, and the minimum ejection 
velocity is of the order of 50 m s”!. These results are fully consistent with 
results of Vokrouhlicky et al. (2006b,c,a), and regarding the velocity-size rela- 
tionship also with Cellino et al. (1999), further strengthening our knowledge 
about the EVF in asteroid families. 

A similar approach is afterwards applied to some other asteroid families. In 
addition to the parameters of the Astrid family described in Section 3, using the 
fact that the vy component of terminal ejection velocities is quite leptokurtic, 
Carruba (2016) has also estimated terminal ejection velocity parameter to be 
the range Vej = 5*1" m s”!. The author excluded values of Vej larger than 
25 m s^! from constraints from the current inclination distribution. Using 
similar arguments and methodology, Carruba et al. (20172) found the terminal 
ejection velocity parameter of the Hoffmeister family. Based on the constraints 
from the observed inclination distribution, the authors concluded that the 
family V g; should be lower than 25 m s^!, with Vej = 2045 m s^! being the 
most likely value. Note that the obtained ejection velocity parameters match 
well the escape velocity from the parent bodies of these families. Additionally, 
independent verification of the size of the initial EVFs was obtained by Bolin 
et al. (2018b), based on the V-shape methodology. 

The list of works on the initial velocity fields of asteroid families given 
here is by no means complete. It just intends to demonstrate the main chal- 
lenges and results obtained along these lines, focusing primarily on the more 
recent papers. Nevertheless, these demonstrate that the past efforts allowed us 
to understand the typical sizes of the initial EVFs and a general size-velocity 
dependence. Especially, the inverse dependency of ejection velocities on frag- 
ment sizes has been verified, as well as the sizes of the initial EVFs, which 
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approximately correspond to the escape velocities from the corresponding 
parent bodies. 

We still do not understand well asymmetries in the EVFs, though available 
evidence suggests they are not uncommon outcomes of the disruption events. 
Just as an illustration, Vokrouhlicky et al. (2006b) found that the velocity 
component normal to the Eos orbit plane is by a factor 4 larger than the 
along-track velocity component. A similar situation was also noted in the Ver- 
itas family (see e.g. Novakovié et al., 2010a). While the theoretical motivation 
for asymmetric velocity fields is well justified (see Bottke et al., 1994), their 
quantification for individual families is challenging. Additional works are also 
needed to understand better how the initial EVFs in asteroid families depend 
on material properties and the internal strength of parent bodies. Future works 
are expected to move the step further, synthesize the results into the global pic- 
ture, and connect them with other relevant aspects such as velocity distribution 
among colliding asteroids or interior structures of asteroids. 


3.2 The Yarkovsky/Y ORP-based age determination 


An underlying fundamental issue for proper understanding of many problems 
that could benefit from the knowledge extracted from families is that apriori 
we do not know the ages of those groups. Dating the asteroid families thus 
comes as a critical topic, and once families are identified, we want to know 
the time of the initial collision that formed them. In principle, ideally, age is 
what we should also know to model the long-term evolution of an asteroid 
family. However, as discussed above, age determination is often coupled with 
modelling efforts. Still, different long-term evolution models present different 
levels of coupling (interdependence) with the age of the family. 

Several dating methods are developed and available, each with its own 
advantages and disadvantages. When dating asteroid families, we need first to 
distinguish between young (and very young) families formed less than 10 Myr 
ago, and those older than about 50 Myr. It is because orbits of members of 
young families could still be accurately tracked backwards until the time of 
the fragmentation event. On the other hand, older families do not offer that 
possibility, and considerably different dating methodologies need to be applied. 
As the young families are extensively discussed in Section 4.1, including their 
appropriate dating methods, this section focuses on the age determination of 
older families. We review here the Yarkovsky/YORP chronology-based dating 
methods, which are the most promising approach for older families. One of 
the method's main advantages is that it can be applied to almost all families 
older than 50 Myr, while providing relatively reliable results. 

How does it work? Asteroid families evolve over time, and this is the key 
to computing the time of the initial collision: the smaller pieces generated by 
the impact tend to change their initial configuration, moving far away from 
the parent body. If we can understand how far and how fast they moved over 
time, we can then find the age of the family. 
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'The Yarkovsky effect is a non-gravitational perturbation produced by the 
anisotropic thermal emission of thermal photons from the asteroid surface. It 
depends on the object's physical properties (e.g. thermal inertia, Bond albedo, 
density) and in particular, it is proportional to the inverse of the diameter. 
As a result, its effect is larger on smaller bodies (Vokrouhlicky et al., 2000; 
Farnocchia et al., 2013), producing a substantial drift of the semi-major axis 
over million of years. Due to the Yarkovsky effect, the smallest members of the 
family tend to move farther away from the parent body, creating the so-called 
V-shape. 

'This behaviour becomes clear when we plot the V-shape of the family in 
the (ap, 1/D) or (ay, H) plane, where a, is the proper semi-major axis, H is 
the absolute magnitude, and D is the diameter of the family members (see 
Fig. 6). 

At the beginning of this century, Bottke et al. (2001) set the stage for the 
Yarkovsky based chronology of asteroid families by recognizing the role of the 
effect in spreading families in terms of the orbital semi-major axis. 

Nesvorny et al. (2003) were the first to apply the Yarkovsky chronology- 
based method to estimate the age of an asteroid family.? In this work Nesvorny 
et al. (2003) used a simple approach by fitting an envelope to the distribution 
of the family in the (ap, H) space. In this way, they were able to estimate the 
age of the Themis family to be 2.5+1.0 Gyr. Later on, Nesvorny et al. (2005b) 
used the same approach to estimate the ages of many other asteroid families. 

A step forward in Yarkovsky chronology-based age estimation was taken 
in papers by Vokrouhlicky et al. (2006b,c). In these works, instead of simply 
fitting an envelope in the (ay, H) plane, the authors used a Monte Carlo- 
like numerical approach. They fit several parameters simultaneously, including 
the initial velocity field of the family and the distribution of family members 
in terms of the semi-major axis. Observation of the coupled influence of the 
Yarkovsky and YORP effects was also a novel aspect of this work. This allowed 
the authors to obtain more robust age estimations for several families. 

A generally similar, but yet somewhat improved, methodology to the one 
used by Nesvorny et al. (2003) was recently proposed by Nesvorny et al. (2015). 
The authors opted to estimate the family ages based on the following equation: 


C a y? p 0.2\? 
Tage 2 1G E AA y ee. G 
s oe ES =) la x) (= >) E ei 


where p is the asteroid bulk density and p, is the visual geometric albedo. The 
constant C is related to the age of the family and depends on the semi-major 
axis drift. It is defined as C = 107 7/5|a, — ap, «|, where ap, is a center of the 
family, and H is an absolute magnitude. 

Another way to approach the problem has been proposed by Spoto et al. 
(2015). It is also based on the finding an envelope to the distribution of the fam- 
ily, but in the (ap, 1/ D) space, and this time the envelope is computed through 


5The idea was for the first time presented a bit earlier by Vokrouhlicky et al. (2002). 
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a least-squares fit. An additional difference concerning described techniques is 
that the age is here estimated from the slopes of the V-shape, not from the 
position of the V-shape borders. We will explain soon why this matters, but 
let us first explain the method in more detail. 

The first step consists of finding how much the orbits of the smaller mem- 
bers moved over time, i.e. how much the proper semi-major axis changed. This 
is obtained through a linear regression fit applied independently to the two 
sides of the V-shape: 


e The two sides of the V-shapes are sampled using bins of different sizes. The 
goalis to obtain bins with a similar density. Consequently, bins will be larger 
at the bottom and smaller at the top of the V-shape. 

* For each bin, a representative is chosen. For the left (inner) side, it corre- 
sponds to the element with the lowest proper semi-major axis value. For the 
right (outer) side, it is the element in the bin with the largest ay. The two 
sets containing all the representatives can be seen as a sample of the inner 
and outer sides of the V. 

* A linear fit is performed on each side, considering a weighting scheme and 
the outlier rejection as in Carpino et al. (2003). See an illustration shown in 
Fig. 6. 


The inverse of the slope obtained as a result of the fit represents the space 
spanned by the family members due to the Yarkovsky effect that we call Aa. 
'This is true for both sides of the V-shape. 

The second step of the method consists of finding the drift rate at which the 
fragments moved over time. Since the Yarkovsky effect has never been directly 
measured from observations for individual objects in the main belt, Spoto 
et al. (2015) rescaled the value detected for the asteroid (101955) Bennu (e.g., 
Chesley et al., 2014; Farnocchia et al., 2021), to obtain reasonable measures for 
each family, according to their known physical and spectral properties (note 
that the corresponding values for Bennu are accurately known). The formula 
used for the rescaling is: 


da EY da v ABennu(1 = Cr ennu) DBennu PBennu cos ^y 1— A 


dt B dt Bennu ya (1 = e?) D Pp COS YBennu UE ÁBennu 


where p is the density chosen for each family, and A the Bond albedo. It is also 
assumed that cos y — +1 for obliquity, extremizing the effect, depending upon 
the inner or outer side of the V-shape. The age of the family is thus given by: 
Tage = Aa/(da/dt). 

Generally speaking, however, the Yarkovsky effect, and therefore the time 
evolution of the V-shapes, may not necessarily scale as 1/D. For instance, 
Delbo’ et al. (2007) suggested that thermal inertia depends on a diameter of 
an asteroid as T = D^?-1. Based on this, Bolin et al. (2018a) proposed that 
Yarkovsky drift scales as (da/dt) c D~%, and consequently a size-dependent 
modification of V-shapes. However, Bolin et al. did not find evidence for this 
hypothesis among 25 investigated main-belt asteroid families. 
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Fig. 6 "V-shape" of the family of (20) Massalia. The figure represents the family in the 
(ap,1/D) plane, where ap is the proper semi-major axis in astronomical units, and D is 
the diameter in km. Blue crosses identify the family members, while black dots are the 
background objects. The two solid black lines are the result of the fit to the sides of the V 
as in (Spoto et al., 2015). 'The blue lines represent the different steps of the fit while trying 
to converge. The red circle around the number (20) gives the position of the parent body, 
while the black circles around blue crosses identify the position of outliers that the fit has 
rejected. The blue stars represent 3/1 resonant objects. 


In addition to the described cases, an interesting aspect of the Yarkovsky 
chronology is met for families whose members are locked inside the power- 
ful mean-motion resonances. The resonances lock the semi-major axis in these 
cases, and the Yarkovsky-induced drift translates into eccentricity. This mech- 
anism produces a V-shape in the (ep, 1/ D) plane, in an analogy to more typical 
V-shapes seen in the (ap, 1/D) plane. By determining appropriate scaling fac- 
tors, the V-shapes in the (ep, 1/D) plane could also be used to determine the 
family ages (Broz and Vokrouhlicky, 2008; Milani et al., 2017). 

Finally, let us recall here that the age of an asteroid family could be approx- 
imately estimated also based solely on the YORP effect. The idea relies on the 
YORP-eye concept introduced by Paolicchi and Knezevié (2016). The YORP- 
eye refers to a depopulated region in the V-shape plot of the family that could 
be diagnostic of the age of the family (Paolicchi et al., 2019). More in par- 
ticular, the position of this area within the V-shape depends on the age of 
the family. Assuming a non-dimensional parameter YORP-age Yage is simply 
given by the equation: Yage = rrA/a?, where 7; is the age of the family (in 
Myr), A is Bond albedo, and a is the semi-major axis (in au). Then, calibrat- 
ing its value on families with known ages, the ages of other families could be 
estimated. We refer readers to the two previous papers for further details on 
the algorithm. The YORP-eye method could be useful to get an independent 
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estimate of families ages. These estimates, however, are usually affected by 
large uncertainties and, often, not unique. 

What are the good and bad sides of different versions of the Yarkovsky-based 
chronology? And what are their general limitations? 

The first thing that should be understood is that the V-shape is gener- 
ally the result of two processes: the initial ejection velocity field and secular 
evolution of the semi-major axis due to the Yarkovsky effect. At first, larger 
fragments tend to be ejected at lower velocities and smaller fragments at higher 
speeds. Thus, the initial configuration after the collision is expected to see the 
smaller fragments located further away from the parent body (Michel et al., 
2015). After that, being more sensitive to the Yarkovsky effect than the larger 
ones, smaller fragments also move even farther away from the parent body. 
Separating the contribution of these two effects is a complex issue. Therefore, 
this, to some degree, affects the accuracy of all the variants of the method. 
The less sensitive and probably the most accurate should, in principle, be the 
method developed by Vokrouhlicky et al. (2006b,c). Its main shortcoming is 
that it cannot be so easily applied to many families, but needs to be used on 
a case-by-case basis. 

Both variants based on an envelope fitting, either in the (ap, HI or 
(ap, 1/D) plane, are affected by the contribution of the initial velocity field, 
which is not taken into account. An important aspect to realize here is that 
method used by Spoto et al. (2015, see also Section 3.1) depends on the ini- 
tial velocity field differently from the classical Yarkovsky chronology-based 
age determination. The point is that Spoto et al. determined ages from the 
slopes of the V-shapes and not from the position of the V-shape borders. The 
slopes do not generally depend on the initial size of the family, which is the 
case in classical V-shape envelope fitting methods. This, however, does not 
mean that the Spoto et al. version of the method is independent of the initial 
velocity field. With the ejection velocities being correlated with the size of the 
fragments, the initial velocity field already produces a V-shape like distribu- 
tion. For this reason, the borders of the V-shape are not vertical lines in the 
(ap, 1/D) plane, and this initial V-shape slope affects the final age estimate. 

Regardless of the method used, the contribution of the initial velocity fields 
to the current distribution of family members decreases with the ages of aster- 
oid families. Therefore it should be less important for families older than 1 
Gyr. Still, the probability of family formation suggests that families formed 
from larger parent bodies are older (Bottke et al., 2005). At the same time, 
ejection velocities are, in principle, also correlated with the size of the par- 
ent bodies. The larger bodies eject on average fragments at larger velocities. 
Therefore, the older families, which are typically also older, are expected to be 
formed in collision events with higher ejection velocities. Therefore, even for 
older families, the initial velocity field could play a role in age determination. 

'The other aspect may also affect the accuracy of the age estimation. The 
Yarkovsky calibration is one of the most limiting aspects. The effect depends 
on many physical properties such as density, surface thermal properties, spin 
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state, etc., unknown for most asteroids. Another important mechanism is the 
YORP effect, which modifies the spin states of asteroids, determining the 
magnitude of the Yarkovsky effect. Finally, the dynamical effects such as close 
encounters with massive asteroids could randomize the semi-major axis to 
some degree, obscuring the V-shape distribution. 

Having all these in mind, we can expect that current family age determi- 
nations are at best accurate on a factor of about 2. This should be considered 
in any analysis based on the estimated family ages. 


3.2.1 Recent Family Age Estimates 


Spoto et al. (2015); Milani et al. (2017) and Nesvorny et al. (2015) provided age 
estimates, or necessary information to directly compute those ages, for more 
than 50 families. The methods used, along with their limitations, are discussed 
in Section 3.2. The main strength of these estimates is that they are obtained 
using uniform methods, which is particularly important in reconstructing a 
collisional history of the main belt. 

'The obtained results are shown in Fig. 7. The ages referring to Nesvorny 
et al. (2015) are computed using Eq. (2) and information given in this work, 
and the same density as in Spoto et al. (2015), to allow more direct comparison. 
There are a few important points to notice here. The results are generally in 
good agreement and compatible from the statistical point of view. In two cases 
where the error bars do not overlap (for families 569, 363) the disagreement is 
caused by different names used for the same family, or the same name for dif- 
ferent families. For instance, (569) Misa contains a younger subfamily, namely 
(15124) 2000 EZ39. While Spoto et al. (2015) under the age of (569) Misa 
provides the result for the older family, Nesvorny et al. (2015) refers to the 
name (569) Misa to the younger subfamily. Since the significant difference in 
the results. Similarly, the difference in age estimates for (363) Padua should 
result from different memberships. However, the exact reason is not apparent, 
and this case might be worth further investigation. 

In any case, most of the results from the two independent estimates are in 
good agreement. This suggests that the ages are relatively reliable. However, 
it should be noted that uncertainties are significant. Let us illustrate this in 
the (24) Themis family. The method developed by Spoto et al. (2015) uses an 
independent fitting of both sides of the V-shape. In the case of the Themis 
family, the authors obtained the age of 2450 + 840 Myr from the inside slope 
and 3780 + 960 Myr from the outside slope. Though formally statistically 
compatible, these two nominal results are not so close each to other. At the 
same time, data from Nesvorny et al. (2015) suggests an age of 38002-1500 Myr. 
Combining these three estimates, we got that the Themis family is 3340 + 
1970 Myr old. It is a much older solution than the usually quoted age of 
2.5 + 1.0 Gyr found by Nesvorny et al. (2003). At the same time, uncertainty 
is huge. Hopefully, a better solution will be found shortly, as Themis is one of 
scientifically the most interesting asteroid families, with a space mission also 
being proposed (Landis et al., 2022). 
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Fig. 7 Ages and their corresponding errorbars for 30 families estimated both by Nesvorny 
et al. (2015) (red circles) and Spoto et al. (2015); Milani et al. (2017) (blue squares). The 
x-axis shows the number of the family namesake asteroid (e.g. 24 for Themis), while the y- 
axis represents the time span of millions of years. Families are sorted from the oldest (top 
left) to the youngest (bottom right). 


Having mentioned two independent age estimates from the in- and out-side 
slopes of the V-shape by Spoto et al. (2015), it is worth mentioning that in some 
cases, this yields different ages for the same dynamical family. Such a situation 
is encountered when the two sides of the V-shape have significantly different 
slopes, likely because the dynamical family results from multiple collisions. 
The best example is the asteroid (4) Vesta. Two different ages have been 
found for the family (Spoto et al., 2015). They could be interpreted as two 
different families, corresponding to the two largest basins observed by the 
Dawn mission, namely Rhea Silvia and Veneneia. The estimated formal ages 
for the two craters are 1 Gyr and 2 Gyr respectively (Marchi et al., 2012), in 
agreement with the two ages found for the family of (4) Vesta by Spoto et al. 
(2015). 

We discussed here only age estimates from Spoto et al. (2015); Milani 
et al. (2017) and Nesvorny et al. (2015), as these works provided the results 
for many families using the uniform method. There are, however, many other 
age estimates, often given for a single or a few families. In many cases, such 
estimates should be of somewhat better accuracy than those obtained semi- 
automatically for a large number of families. It is beyond our scope to mention 
all these works. Let us, however, discuss one striking example that suggests 
even more caution is needed when interpreting family ages. Yang et al. (2020b) 
recently analysed (31) Euphrosyne family and found an age of 280+480 Myr. 
It is significantly younger than the estimate of 1230 + 410 Myr coming from 
Spoto et al. (2015). Previously, the age of the family was also estimated to 
~ 500 Myr, by modelling population of family members interacting with the ve 
secular resonance (Carruba et al., 2014). Though in principle, one may expect 
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those age estimations based on the evolution caused by gravitational interac- 
tions used by Yang et al. (2020b) and Carruba et al. (2014) could somewhat 
underestimate the age, the differences are still immense and raise a warning. 

From a more broad perspective, Fig. 7 also points out the limitation of our 
current knowledge of asteroid families: 


e There is an evident absence of primordial families in the plot. This bias 
appears because of the clustering method used in the family classifica- 
tion (Nesvorny et al., 2015; Deienno et al., 2021). The hierarchical clustering 
methods usually cannot identify very old families, because they are too 
dispersed in the space of orbital elements. 

e The uncertainties are too large to obtain a reliable chronology, and ages 
overlap. These immense uncertainties come from a lack of knowledge of 
asteroids’ physical and spectral properties and from the rescaling used to 
obtain reasonable values of the Yarkovsky effect in the main belt. 


3.2.2 Future age determination perspectives 


We expect significant progress in the precision of family age estimates in the 
future. On one side, discovering new asteroids and improvements in our knowl- 
edge of the physical properties of asteroids will be very useful to obtain more 
accurate ages. On the other hand, existing methods should be further improved 
and extensively tested. The testing should be practical for quantifying their 
uncertainties and identifying the primary sources of the errors. Regarding the 
physical properties, along with data expected from future sky surveys, an 
example of good practice is a step taken by Sergeyev et al. (2021) in mining 
old archival data to find new information. The authors searched for moving 
objects in the SkyMapper Southern Survey and extracted their photometry 
and astrometry from the images. The result extends the known multi-filter 
photometry and taxonomic classification for more than 200 000 objects. 

Recent work by Dell'Oro et al. (2021) is an example of steps that we need 
to follow to understand better the uncertainties and limitations of the family 
age estimation models. These authors investigated the role of physical pro- 
cesses on age estimations. Indeed, tests performed by Dell'Oro et al. (2021) are 
simplistic, focusing only on some of the relevant effects. The approach neglects 
the full effect of dynamical mechanisms or the YORP effect. At the same time, 
though not explicitly stated in the paper, a number of test particles used to 
represent synthetic families seem to be significantly larger than what we have 
in a realistic situation. Combined with neglected effects of resonances, it makes 
V-shape borders fitting unrealistically accurate. Nevertheless, this is a helpful 
first step, which needs to be improved further. 

Regarding the calibration of the Yarkovsky effect, the ESA Gaia mission 
will represent a real game-changer in the estimation of the ages of aster- 
oid families. Originally designed to create the most accurate 3D map of the 
galaxy (Gaia Collaboration et al., 2016, 2018b), launched in December 2013, 
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the Gaia spacecraft also produces highly accurate observations of solar sys- 
tem objects. The second release, in April 2018, showed for the first time the 
improvements that Gaia is going to make to asteroid orbit determination (Gaia 
Collaboration et al., 2018a). The post-fit residuals of Gaia asteroid observa- 
tions reach the mas level for bright asteroids, two orders of magnitude better 
than what we reach from Earth using ground-based observatories. For the 
first time, we can think about having the same accuracy for main-belt aster- 
oids that we obtain for near-Earth asteroids when we observe them with the 
radar. This is crucial for detecting the Yarkovsky effect in the main belt. This 
non-gravitational perturbation can be measured from the astrometry using a 
non-linear least-squares fit, where the effect is one of the parameters to be 
solved-for (Farnocchia et al., 2013). The more accurate the orbit, the higher 
the probability of detecting the Yarkovsky effect on small objects. Gaia finally 
opened this opportunity for the asteroids in the main belt. The goal is to obtain 
at least one measurement for every family, or at least every spectral type. That 
will help solve the problem of large uncertainties in age determination. 


3.3 Families delivering fragments on Earth 


The fact that families play an important role in producing near-Earth aster- 
oids and meteorite flux on the Earth is not new (e.g. Zappalà and Cellino, 
1992; Binzel and Xu, 1993; Morbidelli et al., 1995). As families evolve over 
time, many fragments ecsape from them and are often not anymore recognized 
as members of these families (Novakovié, 2018b). Some of the escaped objects 
are transported towards the near-Earth space. This topic has fascinated sci- 
entists already for about half a century. It is also a natural consequence of 
asteroid family members being a significant fraction of the whole population of 
main-belt objects. Over the years, we have learned a lot about the related pro- 
cesses, but nevertheless, there are still many open questions, which are among 
the central problems in asteroid science. We primarily need to understand: i) 
what is the overall contribution of each family? ii) how the flux from a fam- 
ily evolves with time? and iii) which specific NEAs (meteorites) originate in a 
given family? 

But before answering the above questions, we need to know which fami- 
lies could contribute some members to the NEAs population. In principle, the 
answer is simple, all of them. However, there are two situations to distinguish 
regarding the time evolution of the flux from an individual family. Only fam- 
ilies close to the main transport routes from the main belt to the near-Earth 
region could inject fragments directly onto those routes. The dynamical life- 
times of the objects located inside the strong orbital resonances, identified as 
the primary transport routes of NEAs, are very short, in many cases only a few 
Myr (e.g. Granvik et al., 2018). Therefore, these directly injected fragments 
could be transported quickly, causing many objects to reach the near-Earth 
region in a short time, soon after the formation of the family. Such disruption 
events in the main belt could cause variation in the flux towards the near- 
Earth region. Zappalà et al. (1998) suggested that the formation of the most 
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populous families located close to the transport routes should have produced 
transient episodes of intense craterization of the terrestrial planets. The dura- 
tion of these asteroid showers depends on the resonance involved and on the 
number of injected fragments in each case. On the other hand, families that do 
not inject the fragment directly onto the transport routes send a more stable 
number of objects into the near-Earth region. Therefore, although all fami- 
lies could deliver fragments to the near-Earth region, only those close to the 
powerful resonances are also expected to cause asteroid showers on Earth. 

The motivation of Zappalà et al. (1998) was to explain how short-lived 
asteroid populations inside the powerful resonances are replenished with fresh 
fragments to keep the population of near-Earth objects in a steady state. It 
was, however, immediately apparent that the ejection velocities of the frag- 
ments in family-forming phenomena are too low to inject enough objects into 
resonances. The problem of bringing new objects into the resonances was solved 
soon afterwards by recognizing the importance of the Yarkovsky effect (e.g. 
Farinella and Vokrouhlicky, 1999). Nevertheless, the hypothesis of Zappalà et 
al. about the transient asteroid showers turned out to be correct. 

NEAs population models are all statistical representations of the ‘where- 
abouts’ of NEAs of different sizes. They are built by studying the influx rate of 
NEAs from different sources, estimating via dynamical simulations the mean 
residence time of incoming objects in different parts of near-Earth orbital 
space, and weighing the results against observations (Bottke et al., 2002; 
Granvik et al., 2018). An implicit assumption behind these statistical repre- 
sentations is that the NEA population is close to a steady state, which may 
be true if/when temporal variations of influx are minor. However, the aster- 
oid showers might not fit into this view, and imply that some family forming 
events in the main belt can break assumptions that the NEA population is in 
a steady state. But do we have evidence for such showers? 

One of the ways to constrain the variation in the meteorite flux on Earth, 
is by studying the meteorites that we have in our possessions, as well as dust 
deposits on our planet. Most meteorites that fall today are H and L type ordi- 
nary chondrites. However, the main-belt asteroids best positioned to deliver 
meteorites are LL chondrites (Dunn et al., 2013). By elemental and oxygen- 
isotopic analyses, Heck et al. (2017) studied this contradiction and found 
that the current meteorite flux is dominated by fragments from recent aster- 
oid breakup events and therefore is not representative over longer timescales. 
Additionally, Heck et al. found that the meteorite flux varied over geologi- 
cal time as asteroid disruptions create new fragments that slowly fade away 
from collisional and dynamical evolution. The current flux favours disruption 
events that are larger, younger, and/or highly efficient at delivering material 
to Earth. Similarly, at the other end of the size spectrum, a recent micro- 
meteorites study found changes in the flux composition, possibly due to the 
direct delivery of a meteoroid swarm from the main asteroid belt to Earth 
(Drouard et al., 2019). Additional evidence for impact flux variations comes 
from the work of Mazrouei et al. (2019) who studied the terrestrial impact 
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Fig. 8 An illustration of the main steps in a possible scenario in the dynamical evolution 
of an asteroid family. From a family formation to the delivery of members into near-Earth 
space. The example shows a scenario for the Karma asteroid family after the work by Pavela 
et al. (2021). 


crater record and found that the impact rate increased by a factor of 2.6 about 
290 million years ago. All this evidence suggests that asteroid families play a 
major role in meteorite delivery to the Earth and that their formation may 
cause variation in the Earth’s impact flux. 

Different aspects, including dynamical evolution and transport routes, size 
of the parent body, initial velocity field, and cosmic-ray exposure (CRE) ages, 
need to be considered to establish direct links between the asteroid showers and 
families. Here we will discuss some of the most engaging examples, focusing 
on more recent results. 

One of the first dynamically based links between meteorites and their par- 
ent asteroids was proposed by Nesvorný et al. (2009). The authors identified 
the Gefion asteroid family as a source of the shocked L-type chondrite mete- 
orites, representing two-thirds of all L chondrite falls. To verify this connection, 
Nesvorny et al. followed two different lines of evidence. On one side, they 
searched for possible locations in the main belt capable of transporting mete- 
orites to the Earth on timescales comparable with the most extreme CRE ages 
of the fossil meteorites (<200 kyr; Heck et al., 2004). On the other side, the 
age of marine limestone in southern Sweden, where the fossils of L-type chon- 
drite meteorites and iridium enrichment were found, as well as the ages of five 
large terrestrial craters, are all found to be about 470 Myr. Because of these 
facts, and given the short transport timescale, the authors looked for a group 
of objects formed around 470 Myr ago. Nesvorný et al. found that only the 5:2 
mean motion resonance with Jupiter at 2.823 AU can transport the meteorites 
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over the given short timescale and that the age of a nearby Gefion family plau- 
sibly coincides with the shock event of L chondrites dated 470 Myr. Therefore, 
the authors concluded that the Gefion asteroid family is likely the source of 
the fossil L chondrite meteorites. 

'The perfect example of how an asteroid disruption in the main belt may 
affect the near-Earth environment is the case of the Flora family. Being close 
to the vg secular resonance, the Flora family has been proposed as a significant 
source of a large number of NEOs and LL-type ordinary chondrites meteorites 
(e.g. Binzel et al., 2015, and references therein). What, however, was not well 
constrained were the timing and magnitude of Flora family members striking 
the Earth. Vokrouhlicky et al. (2017b) recently re-accessed the Flora family 
and performed a very detailed study aiming to constrain better the time vari- 
ations in the influx of near-Earth objects from the family. The authors used 
state-of-the-art tools to simulate the long-term collisional and dynamical evo- 
lution of the family and investigated different assumptions about the initial 
dispersion of the fragments. Though there are some uncertainties in the exact 
number of objects and the timing of their arrival in the near-Earth environ- 
ment, Vokrouhlicky et al. (2017b) undoubtedly showed that the flux gradually 
decreased to the present-day situation after an initial spike of objects. There- 
fore, in an interval between about 100-200 Myr after its formation, the Flora 
family sent an increased number of fragments towards the Earth, likely causing 
an asteroid shower-like event. Such a scenario was possible only because many 
fragments produced in the Flora family forming event were injected directly 
into the vg resonance. 

In Section 3 we already discussed the work of Bottke et al. (2015) in light of 
its dynamical simulation advances. However, in this paper, Bottke et al. actu- 
ally presented a systematic analysis of the transport from the five primitive 
inner main-belt asteroid families (Erigone, Sulamitis, Clarissa, Eulalia, and 
New Polana) to the near-Earth region. The authors primarily aimed to identify 
potential source families of asteroid Bennu, a target of NASA's OSIRIS-REx 
sample return mission. By tracking the evolution of millions of test particles 
started near the most likely main-belt source regions for Bennu, Bottke et al. 
(2015) found that Bennu mostly likely arrived into the near-Earth region via 
the vg resonance. The intermediate source Mars-crossing region (IMC) is a 
possible but less likely source region, while the 3:1 resonance is very unlikely 
to transport an asteroid onto a Bennu-like orbit. The applied methodology 
allowed Bottke et al. to estimate that the relative probabilities of Bennu com- 
ing from the vg resonance and the IMC region is 82% and 18%. Using the 
Monte Carlo simulations of the semi-major axis drift, including also for the first 
time stochastic YORP model, it was possible to measure the flux of Bennu- 
like objects from each of five considered families into each relevant transport 
route. Combining the simulations results and the estimates of the probability 
of reaching a Bennu-like orbit via specific transport routes, the authors found 
that the New Polana and Eulalia families had about 7096 and 3096 probabil- 
ity of producing Bennu, respectively. Interestingly, the same methodology was 
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applied to asteroid Ryugu, a target of JAX A's Hayabusa 2 sample return mis- 
sion, and found similar source region probabilities as for a Bennu, opening a 
possibility that both objects originate from the same parent body. 

In more general terms, it is important to underline here that the results 
obtained by Bottke et al. (2015) are size-dependent, meaning that objects 
of different sizes placed on the same orbit as Bennu and Ryugu could have 
different families as the more probable source location. For instance, as noted 
by the authors, for bodies that are considerably smaller than Bennu, the flux 
from families like Erigone should start to play an increasingly important role. 

A link has also been established between aubrite meteorites and the Hun- 
garia family (see e.g. Greenwood et al., 2020, and references therein). However, 
as shown by Cuk et al. (2014), the meteoroids from the Hungaria region pre- 
dominantly reach Earth due to Yarkovsky induced drifting across the orbit of 
Mars, with no assistance from orbital resonances. This results in their longer 
(~ 50 Myr) delivery times, in agreement with the CRE ages of aubrites. There- 
fore, though the Hungaria family definitely supplies meteorites to the Earth, 
it probably did not cause any transient shower of these objects. 

An interesting example is also the Baptistina family. Based on the numeri- 
cal simulations of the collisional and dynamical evolution of this group, Bottke 
et al. (2007) proposed the Baptistina family as a source of Cretaceous/ Ter- 
tiary (K/T) impactor. Despite an extensive simulation effort, at that time, 
the authors were missing a piece of critical information about the composi- 
tion and density of the family members. It prevents them from adequately 
modelling the transport from the family to the Earth. Once the critical data 
about the properties of the Baptistina family became available, it was appar- 
ent that potential Earth's impactors, coming from this family, neither match 
the expected composition of the K/T impactor nor the timing of impacts coin- 
cide with this event (Reddy et al., 2009; Masiero et al., 2012). Nevertheless, 
although the Baptistina is likely not a source of K-T impactor, the work of 
Bottke et al. (2007) still informs us that the family contributes some objects 
to the NEAs population. In particular, it could be an additional source of LL 
chondrite-like NEAs (Vokrouhlicky et al., 2017b). 

There are also other examples of families found to deliver some asteroids in 
the near-Earth region. For instance, Masiero et al. (2015) performed numerical 
simulations of the dynamical evolution of the Euphrosyne asteroid family. The 
family is one of the largest low albedo groups known, and thanks to the large 
orbital inclinations of its members, it is crossed by the vg secular resonance. 
This provides a unique opportunity for an outer main-belt group to deliver 
fragments directly to the vg resonance and via the resonance to the near- 
Earth region. Similarly, Pavela et al. (2021) found that the Karma family, a 
group of primitive asteroids located at the outer edge of the 3:1 MMR with 
Jupiter, which is about 140 Myr old, has been supplying some asteroids to 
the near-Earth region. Currently, there should be about ten family members 
larger than 1 km in diameter, and many more smaller objects, orbiting in near- 
Earth space. finally, Aljbaae et al. (2017) found that the Maria family could 
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contribute some objects to the NEOs population, though the contribution of 
bodies larger than 3 km in size should be practically zero at present. 

The list of families sending objects towards the near-Earth region is cer- 
tainly not complete. We discussed some of the most interesting examples and 
mentioned a few new findings. However, as stated above, any family could, in 
principle, contribute to the NEOs population. Therefore, many new links are 
expected to be established in the future. 


4 Special family classes 


4.1 Young families 


Denomination young in the title of this section refers to a specific type of 
chronology, namely the situation when purely dynamical methods may infer 
the family age. The history of young families started with a seminal paper by 
Milani and Farinella (1994), who studied a compact, outer main-belt Veritas 
family. Analysing the past orbital history of its seven members, they argued 
that this family must have an age of about 50 Myr or less, because two of these 
members -including (490) Veritas- moved away from the tight family zone 
on this timescale. It was the first example when chaotic dynamics were used 
to determine the age of a family. While fascinating by itself, the accuracy of 
the chaotic chronology has its apparent limitations. It would have been more 
promising if a chronology based on regular dynamics was used. A milestone 
in this direction was set by Nesvorny et al. (2002b), who reported a discovery 
of a compact family about the Koronis member (832) Karin, therefore the 
Karin family. More importantly, noting the fairly regular nature of orbital 
evolution in the Koronis zone, they used convergent behaviour of secular angles 
(aka backward integration method - BIM) for 13 orbits of Karin members 
to show that this family is 5.8 + 0.2 Myr old. Little later, Nesvorny et al. 
(2003) re-analysed the structure of the Veritas family in detail, finding that 
this case also offers a possibility of regular-orbit-based determination of its 
origin, aside from the method mentioned above based on chaotic orbits. This 
is because a certain portion of the phase space occupied by Veritas members 
holds sufficiently regular orbits. Nesvorny et al. (2003) thus observed that 
past nodal convergence of Veritas members having regular orbits implies an 
age of 8.3 + 0.5 Myr. The Veritas and Karin families thus became exemplary 
archetypes of young asteroid families. While similar in this sense, they also 
represent quite different worlds, and we comment on them in some detail below. 
The past two decades have seen the astronomy of young families flourishing, 
and we thus overview some of the most interesting results.? 


S'This section presents only some of the most fascinating and best-studied young families. There 
are also other known families younger than about 10 Myr that deserve further study. For instance, 
Novakovié et al. (2012b) present a preliminary analysis of the Lorre cluster and found it was formed 
in a cratering event onto about 30 km large parent body about 1.9 Myr ago. More recently, Carruba 
et al. (2018c) used the backward integration method to study many potentially young families 
and found four of them to be less than 10 Myr old. Tsirvoulis (2019) analysed the population of 
Eos family members and found its young subfamily, namely about 2.9 Myr old Zelima cluster. 
This cluster was later on studied by Carruba and Ribeiro (2020) who used its resonant nature to 
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Karin family.— The discovery of the Karin family by Nesvorny et al. (2002b) 
became an important landmark of planetary science in general. This is because 
aside from recognising a new class of families that are interesting in their own 
right, it helped to solve several standing problems in planetary astronomy. 
Before we get to this context, we first review the main steps in understanding 
the Karin family itself. 

Interestingly, each step in a better definition of the Karin family brought 
along new scientific results. The initial paper of Nesvorny et al. (2002b) had 
39 members, of which 13 were suitable for an experiment on the past conver- 
gence of the secular angles. The propagation model contained gravitational 
perturbations from the planets being the primary driver of orbital precession 
in space. The next step was taken by Nesvorny and Bottke (2004), who made 
use of quickly increasing new asteroid discoveries. They identified already 90 
Karin members, of which 34 had suitable orbits to allow improvement of the 
family age to 5.75 + 0.05 Myr. The smaller uncertainty achieved was not 
the only result of a larger asteroid sample, but it also stemmed from the 
improvement of the dynamical model. Aside from the gravitational perturba- 
tions, it now contained thermal accelerations (the Yarkovsky effect) for each 
Karin member. Determination of the associated secular drift rate of the semi- 
major axis ày was an important side product of the experiment because it 
allowed broadly constraining the surface properties of these small main-belt 
asteroids. T'he possibility to promptly identify family interlopers, otherwise 
a difficult task (see Section 2.2), is an additional important aspect of the 
age-determination via past orbital convergence. This is because the interloper 
orbits do not follow the convergence pattern. As an example, Nesvorny and 
Bottke (2004) proved the originally proposed second-largest Karin member — 
(4507) Petercollins- does not belong to the family. It has massive importance 
in the interpretation and model reconstruction of the family size distribution. 
Carruba et al. (2016b), revisiting the Karin family identification, revealed yet 
another dynamical effect. By that time, they had already 576 Karin members, 
out of which 480 were suitable for repeating the past convergence method of 
Nesvorny and Bottke (2004). The much larger sample of family members for 
which orbits converged had an immediate result of improving the accuracy 
of Karin's age: 5.75 + 0.01 Myr. More importantly, though, Carruba et al. 
(2016b) also estimated and analysed the semi-major drift-rates ày due to the 
Yarkovsky effect for the 480 Karin members with sizes in between 1 and nearly 
6 kilometres. The values of ày for 55 asteroids with sizes in between 2.5 and 
3.5 km peaked at zero value and extended symmetrically to some maximum 
negative and positive values compatible with theoretical expectations. T'hey 
were consistent with a sample of objects having an isotropic distribution of 
spin axes. However, the ày values for 280 asteroids with sizes between 0.9 and 
1.7 km told a different story. Their distribution was also roughly symmetrical, 
but had maxima at some negative and positive values with only a very few 


set independent constraints on the age and initial ejection velocity field. Note also those young 
families connected with main-belt comets are presented in Section 4.3. 
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objects having dy œ 0. Carruba et al. (2016b) showed that this distribution 
is compatible with a sample of objects whose initially isotropic pole distribu- 
tion was tilted by the YORP effect toward ecliptic poles (thus extreme values 
of the obliquity). Therefore the analysis of the Karin family origin brings evi- 
dence of not only the Yarkovsky effect, but also the YORP effect. There is 
another important context for this result. Note that the chronology of mod- 
erately old families proposed by Vokrouhlicky et al. (2006c) makes use of the 
same synergy of the Yarkovsky-YORP effects to reproduce the polarisation of 
small members in the family toward extreme values of the semi-major axis. 
Here the same dynamical processes are caught in action at an earlier phase, 
when the effect on the semi-major axis is still small. However, the bimodal 
obliquity distribution is already apparent. 

Given the pristine state of the Karin family, observed only a few Myr after 
its formation, it became a primary target suitable for comparison with results 
from numerical simulations of collisional fragmentation of large asteroids. Only 
some of the best-characterised very young asteroid families (Sec. 4.2) may rival 
Karin's exclusive status by now. So far, the most detailed work was presented 
by Nesvorny et al. (2006a). These authors found that the Karin family was 
formed by an impact of ~ 5.8 km projectile onto ~ 33 km parent body with 
c 6— 7 km s^! speed and ~ 45° impact angle. The largest created remnant, 
(832) Karin, has about 17 km size, and the size distribution of Karin frag- 
ments then follows from the second largest fragment of ~ 5.5 km with a steep 
cumulative power law of = —5.3 exponent extending to about 2 km sizes. 
Modelling thus indicates that Karin resulted from a rather energetic collision 
with = (0.1 — 0.15) mass ratio between the largest remaining fragment and 
the parent body. As to the ejection velocities, Nesvorny et al. (2006a) found 
that c 3 km fragments (typical to the steep leg in the size distribution) have 
mean barycentric ejection speeds of ~ 12 m s^! with the fastest (and smallest) 
launched at ~ 30 m s-!. These values are comparable, or only slightly larger, 
than the escape speed from the parent asteroid, namely ^ 20 m s”!. While 
fairly complete and satisfactory, Nesvorny et al. (2006a) pointed to several 
unresolved issues that need a more complete study in the future (especially if 
combined with a new update of Karin's population). For instance, predicting 
the right order of magnitude of the ejection velocities of Karin members from 
the largest remnant, their exact 3D velocity field was not successfully repro- 
duced. Related to this issue may also be the fact that the model so far has 
failed to reproduce spin parameters of the largest remnants: (i) (832) Karin is 
a slow prograde rotator with the period of 18.35 hr (e.g., Slivan and Molnar, 
2012), (ii) while the next six to seven large fragments are retrograde rota- 
tors (some of which also have rotation periods determined and found to be 
shorter than that of Karin). Matching such a configuration is still a challenge 
for impact modelling. The spin state of (832) Karin itself may be a sort of 
curiosity, if not a mystery, as it bears similarity with Slivan states of larger 
and Gyr-old members in the Koronis family (e.g., Vokrouhlicky et al., 2003). 
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However, before going too far along with this speculation, one first needs to 
break the longitudinal ambiguity of Karin's pole through future observations. 

While expectedly very useful for the physics of collisional disruptions, 
discovering the Karin and Veritas families, identifying their properties and 
determination of their young age, became crucial in solving the following 
decade-long problem in planetary science. Analysing all-sky scans in ther- 
mal waveband taken by IRAS spacecraft, Low et al. (1984) discovered warm, 
large-scale structures (bands) extending symmetrically up to 10? away from 
the ecliptic plane. The inferred temperatures of 150 — 200 K pointed to their 
origin in the asteroid belt, most likely dust released in asteroid collisions. A 
more complete description of these interesting structures was later presented 
by Sykes (1988), and confirmed by observations of the COBE spacecraft (e.g., 
Reach et al., 1997; Kelsall et al., 1998). Soon after their discovery, sources 
for the dust in these bands were passionately debated. Dermott et al. (1984), 
observing that the latitudes of the most prominent o, 3 and y bands roughly 
coincided with the mean proper inclinations of the prominent Themis, Koro- 
nis and Eos families, proposed that these giant and old families as the sound 
sources. À competing model, in which the dust in the bands is produced in 
Myr-old collisions of moderately large asteroids (15 — 20 km or so), was for- 
mulated by Sykes and Greenberg (1986). While the former hypothesis was 
initially preferred, support for the latter view stemmed from about a degree 
mismatch between the Eos proper inclination and the latitude of the y band, 
and by realising that the age of Themis, Koronis or Eos families is much larger 
than the timescale to reach collisional equilibrium with the whole main belt 
for their dust component. So the traces of their dust should be presumably 
smeared into a smooth signal from the whole belt. Nesvorny et al. (2002b) and 
Nesvorny et al. (2003) immediately recognised the value of the young age of 
the Karin and Veritas families in the context of the discussion on the IRAS 
bands origin. These sources perfectly fit what Sykes and Greenberg (1986) had 
in mind. In fact, the proper inclination of the Veritas family matched very well 
with the @ band latitude, and removed thus the pending problem with Eos. 
A detailed numerical model of the scenario, in which Karin and Veritas were 
the sources of the 8 and y bands, was developed by Nesvorny et al. (2006b). 
Of the many results in this paper, we mention here an important quantita- 
tive justification. Nesvorny et al. (2006b) showed that the observed absolute 
flux from the 8 band could be explained by the total volume of Karin-released 
lum to 1 cm particles equivalent to a sphere of 11 km diameter (and sim- 
ilarly for Veritas and the signal from the y band). This comfortably fits in 
the above-mentioned disruption model of the Karin family parent body. The 
< 10 Myr old Beagle family was little later found to be the source of the a 
band by Nesvorny et al. (2008) (it is possible that the Beagle family is slightly 
older, or holds substructures, e.g., Carruba, 2019, but the connection to the a 
band would not be changed). Strong support for the young-families origin of 
the IRAS dust bands also comes from an independent analysis of Farley et al. 
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(2006), who found traces of the Veritas-released dust in ?He anomaly detected 
in the globe-wide geologic layers of the appropriate age. 

'The second broad-scale topic in planetary science to which discovery of 
Karin (and other young families) significantly contributed has to do with space 
weathering. This phenomenon, known since the analysis of the Apollo-era lunar 
samples, aims —in the asteroid context- to resolve a discrepancy between the 
observed asteroid reflectance spectra and those of the presumed meteorite 
analogues (such as the S-type class of asteroids and the ordinary chondrites). It 
has been suggested that the asteroid reflectance spectra are being modified over 
time due to irradiation by cosmic rays, impingement by micrometeoroids and 
solar wind particles, their combination or other yet to be identified processes 
(see reviews by Clark et al. (2002) or Chapman (2004)). This was known 
by the end of the millennium. However, the characteristic timescale of the 
weathering processes, which may also have the power to discriminate between 
different variants of them, was not constrained by astronomical observations. 
The first attempt in this direction was provided by Jedicke et al. (2004). These 
authors used a compilation of ages for known S-type asteroid families and 
observed their correlation with the average spectral slope from broad-band 
photometry of their members. To make the analysis robust enough, it was 
crucial to compare the data for Gyr-old families (such as Maria, Eunomia 
or Koronis) to young families (such as Karin and Nele/Ianinni, see Nesvorny 
et al. (2003) or Carruba et al. (2018b); later Vernazza et al. (2009) extended 
the sample to even very young families, Sec. 4.2). The extrapolation of the 
colour-age correlation pointed rightly to the ordinary chondrite values even 
beyond the value observed for the young families. This implies very fast onset 
of a weathered surface, while the observed slope of the age-colour correlation 
described the slow tail of the processes ongoing for Gyrs. 

Could we anticipate what comes next? The asteroid database keeps extend- 
ing with discoveries of small asteroids, and this trend may even accelerate in 
the future due to output from powerful sky surveys. Hence, it is perfectly possi- 
ble that the known Karin population could grow to a few thousand in the next 
few years. Observing the trend described above, the Karin age could be readily 
improved. More importantly, though, the completion of the Karin population 
down to a few hundred-meter sizes would allow testing if the size distribution 
keeps increasing along the power-law specified by large members (Nesvorny 
et al., 2006a), or if a slope change occurs. Because a 500 m size asteroid should 
have a collisional lifetime of about 200 Myr (e.g., Bottke et al., 2005), in little 
more than 5 Myr of Karin's age, only about 3% of the fragment population of 
this size should be affected. So if a pronounced slope change is observed above 
this size limit for Karin fragments, it would be more likely related to the frag- 
mentation physics of the parent body. It would be exciting to see what the 
collisional codes predict at this level of resolution. However, there might be an 
unexpected problem to deal with when so many new asteroids are discovered. 
Recall that the primary identification method of the Karin family is based on 
the hierarchical clustering tool in the three-dimensional space of proper orbital 
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elements. With the main-belt population reaching a million asteroids or more, 
the proper element space becomes very crowded. Already the Karin family 
identification, despite its compactness, suffers from overlap with the Koronis-2 
family (a recent cratering event on (158) Koronis; e.g., Molnar and Haegert, 
2009). Separating the two families may become a delicate issue in the future. 
Perhaps the past orbital convergence may be the key tool, but it is yet to be 
tested to what size limit it would provide reliable results. 


Veritas family.— The Veritas family analysis by Milani and Farinella (1994) 
brought a novel point of view on age determination, but it was not very con- 
straining. A follow-up work of Knezevié and Pavlovié (2002) used many more 
orbits of Veritas members and was more detailed in analysis, but in fact, 
resulted in an even looser age constraint of less than 100 Myr. A masterful 
step in using the chaotic chronology was thus presented by Tsiganis et al. 
(2007). These authors first looked in detail at the orbital space occupied by 
the Veritas family, and using Lyapunov times associated with each orbit in this 
zone, Tsiganis et al. identified main zones of interest: (i) two chaos-revealing 
groups A and B, characterised by a short Lyapunov timescale between 10 and 
50 kyr, and (ii) two regions R1 and R2 containing much more stable orbits. 
Importantly, the two largest remnants in the family are located in the oppo- 
site camps: (490) Veritas in size, the most chaotic group A, (1086) Nata in the 
regular region R1 (see already Nesvorny et al., 2003). Focusing mainly on the 
methods of the chaotic chronology, Tsiganis et al. (2007) analysed in detail 
dynamics in zones A and B, which are related to the three-body mean motion 
resonances 5J-25-2 and a higher-order 3J 4-35-2. They observed a spread of the 
group-A and group-B Veritas members in proper eccentricity and inclination 
values. Tsiganis et al. found them larger than would have been expected by 
extrapolation from the adjacent regular zones R1 and R2. The extension was 
naturally interpreted as the chaotic diffusion in the respective resonance over 
time. Thus, by conducting a numerical experiment in which orbits are started 
in the expected initial region inside each of the resonances, one can let the time 
go until the spread of the synthetic orbital population matches the observed 
population of the Veritas members. Using an elegant one-dimensional diffusion 
model, in which the fundamental role is played by the respective action invari- 
ant of the two-body heliocentric orbit, Tsiganis et al. (2007) proved that the 
group-A dynamics needs a timescale of 8.7+1.7 Myr (result that was later even 
improved using the same method to 8.7 4: 1.2 Myr by Novaković et al., 2010a). 
This was in a marvellous agreement with the Veritas age determined from 
dynamics of regular orbits in the region R1 by Nesvorny et al. (2003). Trou- 
ble, though, occurred with the dynamics of the group-B members. Because the 
associated 3J+38-2 resonance is much weaker, an order of magnitude longer 
timescale is needed to explain their spread in proper eccentricity. This prob- 
lem lasts till now. The most likely solution consists of a spurious association of 
background objects to the family by the clustering method in proper element 
space. In fact, this may partly apply to the objects in the group A. Intrigu- 
ingly, (490) Veritas itself was proposed to possibly be such an interloper, as 
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attempts to include it as the largest remnant in numerical reconstructions of 
the Veritas family size distribution did not lead to satisfactory results (e.g., 
Michel et al., 2011). 

As noted above, the possibility to use regular dynamics for Veritas family 
age determination was first discovered by Nesvorny et al. (2003). Tsiganis et al. 
(2007) and Novakovié et al. (2010a) also used this possibility as a comparison 
case to their chaotic chronology efforts. The novelty was that they used now 
also orbits from the R2 region, and Novakovié et al. (2010a) the newly defined 
R3 region (a tail of the Veritas family containing orbits with semi-major axes 
larger than group A members in the 5J-25-2 resonance). However, all these 
works used convergence of the orbital nodes only, deeming the much faster pre- 
cessing perihelia unreliable. A full-frontal approach to using regular dynamics 
for Veritas age determination was finally taken by Carruba et al. (2017b). Out 
of their 705 orbits in the R1, R2 and R3 regular zones, they selected 274 most 
stable ones, and successfully showed their simultaneous convergence in both 
the longitude of node and perihelion. Inspired by the Karin family analysis, 
they used the Yarkovsky semi-major axis drift rates ày as a free parameter for 
each of the integrated orbits. This surplus of free parameters in their numerical 
experiment was the key element of their success. However, they also verified 
that the required ày has the expected order of magnitude. The resulting Veri- 
tas family age was estimated to be Box Myr. Despite these advancements, 
the *complicated world" of the Veritas family keeps holding mysteries (e.g., 
Novakovié, 2012). 


4.2 Very young families 


Very young families are compact clusters of asteroids sharing similar helio- 
centric orbits with ages younger than 1 — 1.5 Myr. Such ages imply that the 
values of secular orbital angles, the longitude of node (2 and perihelion c, still 
have similar values (as well as semi-major axis, eccentricity and inclination). 
'This is because the gradient of their proper frequencies s and g with a semi- 
major axis is typically several tens of arcseconds per year and astronomical 
unit. Given the initial ejection velocities of the order of ten meters per second 
in these structures, the fragments are initially dispersed in a zone of about 
(1—2) x 107? au in the semi-major axis near the largest remnant. So we roughly 
estimate that in a Myr, the secular angles disperse a couple of tens of degrees 
at maximum. Note that this value nearly always entirely dominates the ini- 
tial secular angles dispersal by the ejection velocity field, which is not larger 
than a fraction of a degree (see, e.g., Eq. (1) in Vokrouhlicky et al., 2017a, for 
a more quantitative discussion in the specific case of the Datura family). The 
fact that very young families appear as clusters of asteroids in five-dimensional 
space of osculating or mean orbital elements (excluding mean longitude) con- 
veniently alleviates the problem of phase space filling with an ever-increasing 
number of known asteroids. T'hese structures appear as a cluster in the space 
of proper elements, but with a small number of known members (sometimes 
only three or four), their statistical significance may be poor. Adding two more 
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dimensions represented by the secular angles helps isolate them from the back- 
ground population of asteroids in most cases. Therefore the primary arena in 
which very young families are searched is the above-described five-dimensional 
space of osculating or mean orbital elements. In fact, the usefulness of proper 
orbital elements, in this case, is questionable (or limited at least) because the 
age of these families is often much shorter than the timescale over which the 
synthetic proper elements are constructed. Their values may thus average out 
some family structures of interest. 

While many clustered structures, potentially very young families, have been 
reported in the literature, here we overview only those cases studied in suffi- 
cient detail. It may mean they contain many fragments, were the subject of 
targeted astronomical observations, or are of some particular interest. 


Datura family.— Historically, a compact cluster of six small asteroids about 
(1270) Datura (thence the Datura family) was the first example of a very 
young family discovered (Nesvorny et al., 2006). Its convenient location in 
the innermost part of the main belt, and large orbital eccentricity of = 0.21 
(thence low perihelia), allowed to quickly increase Datura's known population 
to the present-day census of 76 (see gradual steps in Vokrouhlicky et al., 2009; 
Rosaev and Plávalová, 2017; Vokrouhlicky et al., 20172). Fig. 9 illustrates the 
novel aspect of the very young families in the case of Datura, namely the 
confinement of the secular angles to a limited interval of values quite smaller 
than the full range of 360°. The characteristic signature of the very young 
families is also the correlated nature of the secular angle distribution. The 
correlation is linear, as shown in Fig. 9, when a simple shear of secular angles 
due to the semi-major axis dependence of their precession frequencies, and 
their slightly different initial positions in the family, dominates the effect. As a 
result, fragments initially starting with a slightly smaller semi-major axis than 
(1270) Datura are on Fig. 9 driven to the bottom right corner and vice versa. 
Similarly to the young families with regular-enough orbits (such as the Karin 
family), the age-determination of the family is basically a reversal of the same 
process, namely a numerical experiment of backward orbital propagation in 
which the currently observed configuration of secular angles (such as in Fig. 9) 
is shrunk to the expected initial state. Selecting a subset of dynamically stable 
orbits and the largest fragments in the family (the least affected by the thermal 
accelerations), Vokrouhlicky et al. (2009) obtained 530 + 20 kyr for Datura 
(see also Vokrouhlicky et al., 20172). 

'The convenient location of the Datura family, and the interest rose by 
discovering its very young age in 2006, allowed and motivated several dedicated 
astronomical observations. In this sense, it has an exquisite status among the 
groups of very young families to date. 

Early photometric observations of (1270) Datura were obtained by Wis- 
niewski et al. (1997) and Székely et al. (2005), and allowed us to determine 
its short rotation period of 3.35 hr. A dedicated campaign in the late 2000th 
provided enough additional photometric data to resolve the Datura rotation 
pole at = 76? ecliptic latitude and a convex shape model (see Vokrouhlicky 
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Fig. 9 Correlation of osculating values of secular angles, the longitude of node (2 at the 
abscissa and the longitude of perihelion w at the ordinate (at MJD 59400 epoch), for 
currently known population of 76 members in the Datura family: black symbols are the 
57 multi-opposition asteroids (Datura itself shown by a diamond symbol), and 19 single- 
opposition asteroids are shown using grey symbols. Members in the exterior M9/16 resonance 
have the nodal values between ~ 87° and > 93.5? nodal longitude (as indicated by the inter- 
val at the left), members with proper semi-major axis smaller/larger than (1270) Datura 
indicated by arrows on the right). The dashed line approximates the correlation using a lin- 
ear function with a slope of —0.5. The largest dispersal around this trend is seen for objects 
affected by the M9/16 resonance (and some single-opposition objects whose orbits are still 
inaccurate). 
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Fig. 10 Left panel: Size distribution of fragments in the Datura family (the largest rem- 
nant (1270) Datura excluded): (i) the grey symbols are 13 Datura asteroids detected by 
Catalina Sky Survey stations from 2005 to 2012, (ii) dashed lines show two possible solu- 
tions for a debiased population of Datura fragments (assuming a simple broken power law, 
see Vokrouhlicky et al., 2017a), and (iii) the red symbols are the presently known total pop- 
ulation of 76 fragments as of March 2022. Right panel: Shape models for the three largest 
Datura fragments based on optical photometry inversion: (i) (90265) 2003 CL5 on the top 
(labelled 1 on the left panel), (ii) (60151) 1999 UZ6 in the middle (labelled 2 on the left 
panel), and (iii) (89309) 2001 VN36 at the bottom (labelled 3 on the left panel). The left 
two views are perpendicular from the equatorial plane, the right view is pole on. 
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et al., 2009). Many more photometric observations of smaller fragments in 
the Datura family were taken in the past decade. Those concerning the three 
largest of them are interesting (Fig. 10): (90265) 2003 CL5, (60151) 1999 UZ6, 
and (89309) 2001 VN36 have all slow rotational periods of 23.41 hr, 13.88 hr 
and 73.15 hr, respectively, their rotational poles have all ecliptic latitude larger 
than 55°, and their shapes are very elongated (see Vokrouhlicky et al., 20172). 
Because the lightcurve inversion methods allow for reconstructing only a con- 
vex hull of the shapes, some, or all, of them may be contact binaries. On top 
of that, Vokrouhlicky et al. (2017a) used data about detections of 13 Datura 
fragments by Catalina Sky Survey stations from 2005 to 2012, attempting to 
provide a debiased population down to about 200 m in size. While the num- 
ber of objects was still small, and thus the solution could not be strongly 
constrained, they found at least a range of possibilities (assuming a bro- 
ken power-law size distribution). For instance, they predicted a population of 
45 +12 D > 500 m fragments, while currently, as of March 2022, there are 
about 30 of them known. More importantly, the recent years” impressive growth 
of the known Datura population (red symbols at the left panel of Fig. 10) 
seems to exclude some of the shallowest predicted populations (such as the 
one shown using the hashed line at the bottom and left part of Fig. 10). The 
trend seen in the present data may also suggest that the two-slope, power-law 
model of Vokrouhlicky et al. (2017a) might be too simple and should perhaps 
be replaced with more complex parametrizations in future efforts. Overall, 
the information about the Datura family population, both shapes and rota- 
tion states of the largest members, and the size distribution downward the few 
hundred-meter sizes provides a unique challenge for outcomes from numerical 
simulations of asteroid collisions (see, e.g., Henych and Holsapple, 2018, for 
some initial work). 

Spectroscopic observations provided additional information. For instance, 
Mothé-Diniz and Nesvorny (2008) found a progression from S to Sq and Q type 
spectra as we move from (1270) Datura, about 9 km size largest remnant in the 
family, down to kilometre-size small fragments. A more extensive spectroscopic 
study putting the data in the context of the space weathering processes was 
then published by Vernazza et al. (2009). The fact that the space-weathering 
spectral slope indicator was not much different for the Datura family (and a few 
other very young families) and the young Karin family implies a rapid timescale 
for space weathering onset (possibly related to solar wind implantation in the 
surface matrix). 

Datura family's very young age may qualify it well as a possible source of 
observable interplanetary dust (the same as much larger Karin, Veritas and 
Beagle families; Sec.4.1). Vokrouhlicky et al. (2008) attempted to model both 
orbital and collisional dynamics of dust released during the Datura family for- 
mation. Indeed, they found that Datura may potentially be a source of C and 
D dust trails observed by the IRAS spacecraft (e.g., Sykes, 1988). Further con- 
firmation using both modelling and independent observations will be needed in 
any case. In fact, Nesvorny et al. (2006) imagined more spectacular evidence of 
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the very-young-families dust, namely its direct detection in the Antarctic and 
Greenland ice core drilling (after it accreted onto the Earth). Later estimates 
for Datura and other known families, however, indicated that the amount of 
possibly accreted dust is much smaller than would have been needed. 

Finally, despite the compactness of the Datura family, the orbital his- 
tory of some of its members happens to interact with a weak, high-order, 
mean motion resonance M9/16 with Mars (see already Nesvorny et al., 2006; 
Nesvorny and Vokrouhlicky, 2006). Various aspects and illustrations of the res- 
onant dynamics in M9/16 has been studied by Rosaev and Plávalová (2017) 
and more recently by Plávalová and Rosaev (2021). In general, such orbits 
present dynamical chaos effects that undermine their use for the family age 
determination based on regular-orbits convergence. The resonances crossing 
the very young families, such as the Datura family, are, however, too weak to 
allow chaotic chronology for independent verification of their age (such as in 
the Veritas family; Sec. 4.1). 


Schulhof family.- Another aspect that occasionally could turn a very young 
family analysis into a detective story, can be illustrated in the case of the 
Schulhof family. Pravec and Vokrouhlicky (2009) discovered a small cluster of 
four asteroids about the largest body (81337) 2000 GP37 as a side product of 
their systematic search for asteroid pairs. Soon afterwards, Vokrouhlicky and 
Nesvorny (2011) realized that this is only a sub-cluster of a larger structure 
related to the largest remnant (2384) Schulhof, counting eight members alto- 
gether. (81337) 2000 GP37 with its group resides slightly offset from (2384) 
Schulhof in the orbital space. Such a situation perhaps requires special ini- 
tial ejection velocity field properties, possibly related to the geometry of the 
initial, family-forming impact onto the parent object. Nevertheless, their com- 
mon orbital convergence to (2384) Schulhof some 780 + 100 kyr ago justified 
a union of the two clusters of asteroids. However, Vokrouhlicky and Nesvorny 
(2011) opened Pandora's box by noting that two asteroids, tightly accompany- 
ing (2384) Schulhof in its newly discovered cluster, tend to orbitally converge 
less than 100 kyr. They speculated that the first, larger-scale family-forming 
event had left (2384) Schulhof in an unstable state, facilitating the formation 
of a later secondary sub-cluster. 

While later resurrected by Fatka et al. (2020), and perhaps best exemplified 
in the case of the Emilkowalski family, the idea was retracted by a more thor- 
ough work of Vokrouhlicky et al. (2016). These authors found yet another four 
new members in the family (the current census as of November 2021 is 26), but 
more importantly, they proved that the previously suggested young ages for 
some members were caused by their inaccurate orbits. The take-away message 
is that one has to carefully choose suitable members used for past orbital- 
convergence experiments, avoiding not only the chaos-affected resonant cases, 
but also those residing on single-opposition, or otherwise loosely constrained, 
orbits. 


Adelaide family.- Novakovié and Radovié (2019) reported a discovery of a 
new, very compact cluster about the inner main-belt asteroid (525) Adelaide. 


Springer Nature 2021 IATEX template 


Article Title 51 
D (km) D (km) 
10 5 2 1 0.5 Š 2 1 95 
EN i sp T T T 
g | Adelaide Hobson 
o 
N 
e 
N 
= e OL 4 
ToL al I = 
v= sA 
z 2 ag 
Dei 
S eg 
SS 2 een MT. 
E s / 
Et Fe 1 1 fi fi fi L d s rte 18777 Í fi L d 
12 13 14 15 16 17 18 19 20 15 16 17 18 19 
H (mag) H (mag) 


Fig. 11 Left panel: Cumulative absolute magnitude distribution of the identified 72 Ade- 
laide family members. The largest remnant, (525) Adelaide, is shown by a diamond symbol, 
and small fragments by circles. The upper abscissa converts the absolute magnitudes H to 
size D using geometric albedo py = 0.22, determined for (525) Adelaide by observations 
provided by WISE spacecraft. The grey line, approximating the small-fragment distribu- 
tion, is a power-law N(< H) x 107% with y = 2. This is among the steepest progressions 
ever observed in collisionally-born families. Right panel: Cumulative absolute magnitude 
distribution of the identified 51 Hobson family members. The largest objects, (18777) Hob- 
son and (57738) 2001 UZ160 are highlighted by red and blue diamonds, small fragments by 
circles. The upper abscissa converts the absolute magnitudes H to size D using geometric 
albedo py = 0.2, appropriate to the S-type classification of both largest asteroids. The grey 
line,approximating the small-fragment distribution, is a power-law N(< H) œ 107% with 
y = 0.85. 


Observing convergence of nodal longitudes, they suggested an age of ~ 500 kyr. 
More recently, Vokrouhlicky et al. (2021a) revisited the case of the Adelaide 
family for two reasons: (i) its hinted age appeared similar to that of the Datura 
family, and (ii) the two families are very close to each other in space of proper 
elements (consider the values ap = 2.2347 au, ey = 0.1535 and sin Jp = 0.0920 
for (1270) Datura and ap = 2.2452 au, ep = 0.1487 and sin Ip = 0.1170 for 
(525) Adelaide). Vokrouhlicky et al. (2021a) thus investigated a possibility of 
a genetic relation between the two families, such that fragments created in 
the formation-event of one would trigger the formation of the other. First, 
considering a subset of Adelaide members residing on suitable orbits, they 
confirmed the Adelaide family age of 536 + 12 kyr. However, Vokrouhlicky 
et al. found the case of a causal relationship between the formation events 
for the Datura and Adelaide families unjustified. In all likelihood, the two 
families were formed by an impact of a background-population projectile and, 
by chance, happen to be close to each other in time and orbital space. 
Nevertheless, the analysis of Vokrouhlicky et al. (2021a) pointed out yet 
another interesting aspect. The left panel in Fig. 11 shows the cumula- 
tive magnitude distribution of the Adelaide family members: (525) Adelaide 
appears to be about 10 km size largest remnant in the family, accompanied 
by a large number of sub-kilometre fragments. While it is still possible that 
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some intermediate-size fragments are missing due to observational biases, the 
chances are not too high (see, e.g., population completeness study of Hendler 
and Malhotra, 2020). Approximating the small-fragments tail using a power- 
law N(« H) « 107%, we have y ~ 2. Assuming constant albedo value, this 
translates to a cumulative size distribution N(> D) x D^, with 8 ~ —10. 
Either of the coefficients witnesses an extremely steeply-raising population 
(compare, e.g., with relevant simulations of Ševeček et al., 2019). Undoubt- 
edly, the Adelaide family results from a huge cratering event on (525) Adelaide 
itself. However, the size-distribution steepness of the fragments is quite larger 
than observed among cratering events of large main-belt bodies (e.g., Masiero 
et al., 2013b). Hence again, a good example of motivation for impact-cratering 
modellers to test the codes, and relevant parameter space, in the regime of not 
too strong gravity. An interesting implication of the small-fragment population 
steepness is that the known population of the Adelaide family is among the 
most rapidly increasing over the past period. For instance, data published by 
Vokrouhlicky et al. (2021a) had 52 recognized members as of February 2021, 
while the present-day count is 72 (March 2022; Fig. 11). 


Hobson family.- Pravec and Vokrouhlicky (2009) reported a small cluster of 
four asteroids about the middle main-belt object (57738) 2001 UZ160, possibly 
a very young family with an age of less than 500 kyr. T'hey also noted a nearby 
asteroid (18777) Hobson, but they were puzzled about its association with the 
cluster because of its similar size to (57738) 2001 UZ160. Later on Rosaev 
and Plávalová (2016, 2017, 2018) in a series of papers clarified this issue by 
proving that (18777) Hobson orbit converges very well to the members of the 
2001 UZ160 cluster and must be considered a part of it. Because of its slightly 
larger size, the cluster was renamed the Hobson family. Additionally, by 2018 
Rosaev and Plávalová updated the membership to nine members, noted a 
slight chaoticity of some orbits in this cluster (including that of 18777 Hobson, 
possibly due to the influence of a weak, three-body resonance 9J-8$-2), and 
used a simple past convergence method to infer an age of 365 + 67 kyr. This 
age was independently confirmed by Pravec et al. (2018), who also updated 
the Hobson family population to 11 objects. 

The puzzling issue of a similar-size largest remnants in the Hobson fam- 
ily was reanalysed in a new light by Vokrouhlicky et al. (2021b). First, these 
authors identified asteroids associated with this family from the most updated 
orbital catalogue. Vokrouhlicky et al. found the population increased dramat- 
ically (they had 45 members as of July 2021, a number that even increased 
to the total of 51 by the current date). They noted the extreme compact- 
ness of the Hobson family in the osculating values of their secular angles, with 
only ~ 1? dispersion in longitude of node and ~ 2? dispersion in longitude of 
perihelion (compare with data for the Datura family shown in Fig. 9). There- 
fore, the family membership is statistically robust, and all identified asteroids 
exhibit excellent orbital convergence among each other some 330 kyr ago. The 
cumulative absolute magnitude distribution N(« H) of the resulting sample is 
shown in the right panel of Fig. 11. The conjoint existence of nearly equal-size 
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Fig. 12 Left panel: Illustration of convergence of the mean values of secular angles, the 
longitude of node Q in red and longitude of perihelion w in blue, for a selection of 25 Rampo 
family members. The time at abscissa goes to the past, and the ordinate shows the difference 
of the respective angle with respect to the orbit of (10321) Rampo. Mutual convergence at 
~ 800 kyr is apparent. Right panel: Cumulative absolute magnitude distribution of the iden- 
tified 37 Rampo family members. The largest remnant, (10321) Rampo shown by a diamond 
symbol, small fragments by circles. The upper abscissa converts the absolute magnitudes 
H to size D using geometric albedo py = 0.2, appropriate for the inferred S-type spectral 
classification. The grey line, approximating the small-fragment distribution, is a power-law 
N(« H) x 107% with y = 1.15. 


largest remnants, followed by a disjoint population of smaller fragments, is not 
seen in the case of traditional asteroid families (or can be easily excluded as 
an interloper situation, e.g., Masiero et al., 2013b). It has also not been pre- 
dicted by computer simulations of asteroid collisions so far (see Ševeček et al., 
2017, 2019, for relevant parent-body size). 

By extending the efforts of Sevecek et al. (2019), Vokrouhlicky et al. (2021b) 
in fact, found a collisional regime that allows reproducing the observed size 
distribution of asteroids in the Hobson family. Nevertheless, the impact param- 
eters must be fine-tuned, which decreases the statistical weight of such a 
solution. However, there is yet another interesting possibility. Observing results 
from Pravec et al. (2016), who characterised the abundant population of bina- 
ries among € 15 km size main-belt asteroids, Vokrouhlicky et al. (2021b) 
suggested a scenario in which the parent body of the Hobson family would 
be one of these binary systems. In this model, either the primary or the sec- 
ondary undergoes collisional disruption, leaving the second component intact. 
'The largest collisional remnant, and the former intact component in the binary, 
are thus the two large objects in the Hobson family. Because some 15 + 4% of 
small main-belt asteroids are binaries (Pravec et al., 2016), we statistically see 
more such configurations when a larger sample of very young asteroid families 
are discovered in the future. 

Excellent convergence of Hobson-member orbits in the past also provides 
an interesting constraint of the characteristic dispersal velocity of the observed 
fragments at origin, namely € 10 m/s. This is a very low value, comparable 
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with the escape velocity from the proposed systems in both above-mentioned 
models, and well-reproduced in their simulations. 


Rampo family.— In order to demonstrate that quite more is at hand already at 
this moment, we give short information about the Rampo cluster. Its core, con- 
sisting of two small asteroids near the inner main-belt object (10321) Rampo, 
has been discovered by Pravec and Vokrouhlicky (2009). Being interested in 
asteroid pairs rather than clusters, these authors did not pay much attention 
to it other than noting a probable age between 0.5 and 1.1 Myr. More recently, 
Pravec et al. (2018) found another four Rampo members and improved the 
age using past mutual convergence of their secular angles to 7804430 kyr. 

Here we revisited the status of the Rampo cluster by searching in the 

up-to-date catalogue MPCORB.DAT of asteroid orbits provided by the Minor 
Planet Center (we used the version as of March 2022). We found 37 aster- 
oids, including (10321) Rampo, associated with the Rampo family. Only one 
orbit is single-opposition, but some of the multi-opposition orbits also have 
only a few observations, and they are not very accurate (though enough for 
the association with the family). This represents a great population leap com- 
pared to status in 2018, indicating that the past few years had a significant 
increase in known asteroid population propagated to the population status 
of several very young families. The left panel of Fig. 12 illustrates the past 
convergence of mean values of secular angles for a selection of 25 accurately 
determined orbits. We only used nominal orbits and a simple model in which 
only gravitational perturbations from planets, Ceres, Vesta and Pallas were 
included. All orbits tend to converge some 800 kyr ago, in good agreement 
with results in Pravec et al. (2018). However, more quantitative results would 
need to use more intense numerical effort by considering (i) clone variants of 
the propagated orbits (reflecting initial orbital uncertainty), and (ii) thermal 
accelerations in the propagation model. The right panel in Fig. 12 shows the 
absolute magnitude distribution of Rampo-family members. It appears to be 
well compatible with simulations of the collisional breakup of a single parent 
body, in this case, (e.g., Ševeček et al., 2017, 2019). 
More on the dashboard.— The previous examples of very young asteroid families 
showed situations in which the population of known members swung from 
few to many in the past few years. This population increase typically brings 
along the discovery of some unexpected and interesting features. Since this 
trend is expected to continue, perhaps even accelerate, in the coming years 
due to observations of powerful sky surveys, a lot of exciting asteroid science 
may be expected in this field. We finish this section by commenting on three 
population-scale recent results or trends. 

First, Pravec et al. (2018) brought up an interesting conjecture about the 
formation process of very young asteroid families. The traditional asteroid 
families are exclusively formed by collisional disruption of a parent object, 
once existing on a similar heliocentric orbit in the main belt. Realizing that 
the characteristic size of the largest remnants in the known very young aster- 
oid families is D < 15 km, and observing that it overlaps with the size of the 
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parent body of asteroid pairs (e.g., Vokrouhlicky and Nesvorny, 2008; Pravec 
and Vokrouhlicky, 2009; Pravec et al., 2010, 2019), Pravec et al. (2018) inves- 
tigated whether their formation processes may also be the same. Since the 
asteroid pairs have been convincingly shown to form by rotational fission of 
the parent body prevalently (e.g., Pravec et al., 2010, 2019), the question is 
whether all, or at least some, very young families form the same way. Using 
the same approach, namely correlation between the rotational period of the 
largest remnant in the family versus the cumulative mass of smaller fragments, 
Pravec et al. (2018) argued that the very young families share many similar- 
ities with the pairs. This would support their common formation mechanism. 
While very interesting, this result needs further investigation. In some cases, 
small-fragment populations might have been underestimated in this study due 
to their observational biases. A clear example is the above-mentioned exam- 
ple of the Rampo family. In other situations, the rotational-fission formation 
may hold very well. Consider, for instance, the small cluster consisting of two 
small asteroids about (21509) Lucascavin discovered already by Nesvorny and 
Vokrouhlicky (2006) which still has the same population today. Either no more 
fragments exist in this cluster, or they are so small that their contribution to 
the total fragment size will be negligible. 

Second, Fatka et al. (2020) revisited the issue of multiple formation events 
in the very young asteroid families (postulated, but then retracted, in the case 
of the Schulhof family discussed above). Indeed, they found four candidate 
examples in which possibly two events in recent history might have launched 
different sets of fragments associated with the family. The most interesting is 
the case of the Emilkowalski family. Nesvorny and Vokrouhlicky (2006) discov- 
ered that the middle main-belt asteroid (14627) Emilkowalski is accompanied 
by two small fragments forming thus a young cluster. Convergence of secular 
angles hinted at the age of 220 + 30 kyr, although the orbit of one of the two 
fragments —(126761) 2002 DW10- indicated a mismatch in the convergence. 
The very young age of ~ 200 kyr is, however, a perfect match in a broader 
context set by analysis of the solar system dust bands discovered by the IRAS 
spacecraft (e.g. Sykes, 1988). Indeed, a careful analysis of the 17° M/N bands 
by Espy et al. (2009) and Espy Kehoe et al. (2015) indicated they are incom- 
plete, caught in the process of full formation (see, also, Vokrouhlicky et al., 
2008). This constrains the possible age to less than ~ 300 kyr, and the proper 
orbital inclination 17.2° of (14627) Emilkowalski is just needed. However, a 
more complete analysis of Fatka et al. (2020) revealed a more complex struc- 
ture: (i) two of the confirmed Emilkowalski members —(224559) 2005 WU178 
and (256124) 2006 UK337- indeed converge to the largest remnant in the 
family some 250 kyr ago, but (ii) other four fragments —including (126761) 
2002 DW10- tend to converge approximately 1500 kyr ago (with large uncer- 
tainty, though). So were there two cluster-formation events associated with 
(14627) Emilkowalski, or are some incorrectly modelled effects (including pos- 
sibly the orbital uncertainties) fooling the results? This interesting problem 
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would need to be revisited when more family members are found. A broader- 
context hypothesis has been recently published by Carruba et al. (2020c), who 
argued that asteroid clusters that originated in rotational fission of a parent 
body are more frequently found in young families rather than old ones. While 
interesting, this conjecture also awaits further, more detailed analysis. Note, 
for instance, that the most logical test case -the Karin family (Sec. 4.1)- does 
not seem to support these conclusions (see Carruba et al., 2020b). 

Third, we mentioned above that the known asteroid population had seen a 
great increase over the past few years. As a result, some compact clusters with 
few members now have tens of them. The same reasoning may mean that some 
of the previously found asteroid pairs (e.g., Vokrouhlicky and Nesvorny, 2008; 
Milani et al., 2010; Pravec et al., 2019), configurations of just two asteroids 
sharing very similar heliocentric orbits, may soon move to the category of 
clusters (in fact Vokrouhlicky and Nesvorny, 2008, expected that some of the 
discovered pairs might be just a tip of the iceberg, hiding an asteroid cluster 
below the observational limit). Examples could be already seen, but a more 
systematic search awaits for the future. For instance, Pravec et al. (2019) noted 
that the asteroid pair (4765) Wasserburg and (350716) 2001 XO105 could be 
accompanied by a third body 2016 GL253. In fact, the present-day catalogue 
of asteroid orbits reveals that the Wasserburg cluster contains already six 
members, all showing very good past orbital convergence. In the same way, the 
young asteroid pair (5026) Martes and 2005 WW113 has a third companion 
in 2010 TB155. Altogether, they form a three-object Martes cluster, which is 
presently the youngest in its category (having an age of 18 + 1 kyr; Pravec 
et al., 2019). 


4.3 Water Content in Asteroid Families 


Our view about the water-ice content possibly still present in the main aster- 
oid belt has changed dramatically over the last decade or so. The asteroids are 
mostly considered dry objects, in contrast to the comets, which are expected 
to contain a significant amount of water-ice. This view started to change when 
Hsieh and Jewitt (2006) presented data showing the existence of a popu- 
lation of comets in the main asteroid belt. T'hese objects, coined main-belt 
comets (MBCs), have both the orbital characteristics of asteroids and the 
physical characteristics of comets. MBCs are today a subpopulation of the 
so-called active asteroids, which are asteroids that at least occasionally dis- 
play a comet-like appearance. Different possible mechanisms are initiating such 
manifestation of the active asteroids, such as dust ejection due to impacts or 
rotational fission. However, the activity of MBCs is driven by the sublimation 
of volatile ices, implying that these objects contain some amount of water-ice. 

Spectroscopic evidence of hydration is also common among primitive aster- 
oids, with the most direct evidence being the presence of a strong absorption 
feature around the 3 um band. This feature corresponds to different modes of 
OH groups in hydrated minerals and H20 molecules on the surface of the small 
bodies (Rivkin et al., 2015). Based on the shape of the 3 um band, the asteroids 
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exhibiting this feature could be divided into several categories corresponding to 
the prototype body. The first category (sharp or Pallas-like) exhibits a sharp 
3 um feature, attributed to hydrated minerals. The second group ( Ceres-like) 
exhibits a 3 wm feature with a band centre at around 3.05 um, superimposed 
on a broader absorption feature from 2.8 to 3.7 wm. The third group ( Europa- 
like) exhibits a 3 jum feature with a band centre at 3.15 um. The fourth group 
(rounded or Themis-like) is characterized by a rounded 3 um band (Takir and 
Emery, 2012; Rivkin et al., 2015). All these are, to some degree, hydrated 
objects, with the fourth category being explicitly linked to the presence of H2O 
ice. The representative members of the rounded 3 um feature are asteroids 
(24) Themis and (65) Cybele (Takir and Emery, 2012), with a similar band 
shape also found in members of Themis, Hygiea, and Ursula asteroid families 
(De Prá et al., 2020a). 

Where are the families in this story? The asteroid families are generally 
homogeneous in composition, and therefore the chemical properties of members 
of an asteroid family are expected to be similar. It implies that if one member 
contains water-ice, this is likely the case for the other members. Interestingly, 
the first asteroid with surface ice detected, namely (24) Themis, is a member 
of an asteroid family. Similarly, available spectra of the family members show 
that most of them have a rounded 3 um feature indicative of water-ice. Most 
of the known MBCs also belong to families. 

Therefore, the families are deeply relevant not only for understanding the 
origin of MBCs, but also for the origin and evolution of the water-ice in the 
main asteroid belt. As soon as it was recognized that MBCs typically belong 
to the families, links between the two populations started to be investigated. 
A possible connection between the MBCs and young asteroid families was first 
proposed by Nesvorny et al. (2008). The motivation to link MBCs to young 
families was based on the expectation that ice could survive at heliocentric 
distances « 3.3 au over the age of the solar system, only if deeply buried inside 
the larger asteroids. In this respect, the formation of young families provided 
a natural explanation of how the deeply buried ice has been brought close to 
the surface. In years following the work by Nesvorny et al. (2008), a couple of 
other connections between the young families and MBCs have been proposed 
(e.g., Novakovié et al., 2012a, 2014). 

As the available data has grown over the years, Hsieh et al. (2018) per- 
formed the first comprehensive analysis of the links between MBCs and 
asteroid families. This work found a strong indication that MBCs are primar- 
ily members of primitive asteroid families, regardless of how old these families 
are. Such results suggest, on one side, that parent bodies of primitive families 
were typically water-bearing asteroids. On the other side, these findings imply 
that ice could survive in main-belt asteroids over the age of the solar system 
under certain conditions, in agreement with estimates of the buried ice loss 
rate (Schórghofer, 2008; Schórghofer et al., 2020). Therefore, the young age of 
families hosting MBCs may not be critical, as initially suggested. Nevertheless, 
primitive young families are still plausible sources of MBCs. The more recent 
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associations of MBCs to families further support this view. Novakovié (2018a) 
found that 427P/ ATLAS (2017 S5) is a member of the Theobalda family, while 
P/2017 S8 (PANSTARRS) is the first MBC associated with the Pallas fam- 
ily, and P/2019 A7 (PANSTARRS) is the first MBC known to belong to the 
Luthera family (Novakovié and Hsieh, 2019). Therefore, the evidence support- 
ing the link between the MBCs and primitive asteroid families is growing, and 
it is now very robust.” 

While we still do not understand the links between the families and MBCs 
properly, three families stand out as the main reservoirs of MBCs. These are 
'Themis, Lixiaohua and Theobalda families. These three families are of very 
different ages, but their members all share primitive composition characteris- 
tics. An important step forward in fully solving the mystery of MBCs will be a 
better characterisation of these three families, particularly the Theobalda fam- 
ily, which has been poorly studied so far. T'his is a key to understanding why 
these families are reservoirs of MBCs, while some other primitive families are 
not. In what follows, we review the main facts about the Themis, Lixiaohua 
and Theobalda families and discuss some other families expected to represent 
a birthplace for some MBCs, though still not recognised as such. 


4.3.1 Main properties of three families associated with the 
most MBCs 


Themis family.- The family is one of the first five identified, the so-called 
Hirayama families (Hirayama, 1918). It is a prominent group of asteroids 
located in the outer main belt, with low inclinations and low eccentricities. 
'The Themis family is dominated by compositionally primitive, C-type aster- 
oids, with many of them exhibiting spectra suggestive of aqueously altered 
mineralogy, similar to carbonaceous chondrite meteorites (Florczak et al., 
1999; Landsman et al., 2016). The family mean geometric albedo based on 
the data from Mainzer et al. (2016) is 0.068, while recently Jiang and Ji 
(2021) found that the average thermal inertia of family members is about 
40 J m~? s712 K=, 

The Themis family parent body is believed to be differentiated icy-asteroid, 
about 400 km in diameter (Castillo-Rogez and Schmidt, 2010; Marsset 
et al., 2016), which was catastrophically disrupted roughly about 2.5-3.5 Gyr 
(Nesvorný et al., 2003; Spoto et al., 2015, see also Section 3.2.1). Observa- 
tions and modelling have implied that asteroid (24) Themis is the core of an 
icy planetesimal, and therefore the family members are likely also icy objects. 
The evidence supporting this hypothesis includes a near-infrared absorption 
feature attributed to water ice frost was detected in spectra of both, aster- 
oid (24) Themis (Campins et al., 2010; Rivkin and Emery, 2010) and another 
prominent member of the family, (90) Antiope (Hargrove et al., 2015). 

The population of active asteroids includes at least three Themis 
family members: 133P/Elst-Pizarro, 176P/LINEAR, and 288P/(300163) 


TFor more general links between active asteroids and asteroid families, we refer readers to Hsieh 
et al. (2018). 
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2006 VWi39. A fourth active asteroid, 238P /Read, is considered a possible for- 
mer Themis family member whose orbit has dynamically evolved to the point 
at which it cannot be directly associated with the family (Haghighipour, 2009). 
All these four objects are likely MBCs based on the results of their dust mod- 
ellings and recurrent activities (see Hsieh et al., 2018, and references therein). 
Additionally, the proper orbital elements of a newly discovered active asteroid 
(248370) 2005 QNi73 (Hsieh et al., 2021), which is currently inside the 11:5 
MMR with Jupiter, are also compatible with the Themis family membership, 
though a cause of its activity is yet to be confirmed. 

'The Themis family also includes two known younger subfamilies, the Beagle 
and 288P families. Both subfamilies have been dated within the past 10 Myr 
(Nesvorny et al., 2008; Novakovié et al., 2012a), but more recently Carruba 
(2019) found that the Beagle family could be somewhat older (about 35 Myr). 
In addition to the Themis family, two MBCs have explicitly been linked to 
young families, 133P/Elst-Pizzaro to the Beagle group and 288P to the small 
cluster named after it. 

While there are many ice-rich objects beyond the orbit of Jupiter, the 
Themis family is likely the largest reservoir of icy objects in the asteroid belt. 
As such, it is an important keystone for understanding how icy asteroids form 
and evolve. For that reasons, Landis et al. (2021) recommend the asteroid (24) 
'Themis and its family members as prime targets for further study and suggest 
both telescopic and especially spacecraft campaigns be used to explore these 
bodies in the next decade. 

Because the Themis family is next to the 2:1 MMR with Jupiter, it is 
believed that some family members may have been captured and scattered 
by the resonance since the family formation (Morbidelli et al., 1995). Recent 
work by Hsieh et al. (2020) suggests that there could be Themis family objects 
evolving onto Jupiter Family Comets-like orbits in the present day. Similarly, 
a fraction of family members has potentially reached the near-Earth region, 
and recent work by Schórghofer et al. (2020) suggests that some amount of 
water-ice could survive such transport, representing possibly water-bearing 
near-Earth objects. 


Lixiaohua family.- Like the Themis, the Lixiaohua is also a primitive asteroid 
family located in the outer main belt. The Lixiaohua is, however, somewhat 
smaller and significantly younger than the Themis family. The Lixiaohua fam- 
ily was formed by a super-catastrophic disruption of up to 200 km large parent 
body (Durda et al., 2007), that had happen 1554-36 Myr ago (Novaković et al., 
20102). 

The family mean geometric albedo is p, — 0.044 (Masiero et al., 2015), 
which is about 296 smaller in absolute terms than in the case of the Themis 
family. The spectroscopic properties of the Lixiaohua family are recently stud- 
ied by De Prá et al. (2020b). They found that a majority of the objects belong 
to primitive classes, such as the C- and X-complex, and T- and D-types. In 
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particular, the two largest members, asteroids (3330) Gantrisch and (3556) Lix- 
iaohua, are classified as T-type objects, which is an intermediate class between 
X and D-types. 

De Prá et al. (2020b) also found that Lixiaohua has distinct properties 
from the Themis family (see also Section 4.3.2), which can be interpreted as 
a difference in composition or differentiation level. An important dissimilarity 
between these two groups is related to the 0.7 wm hydration band. While the 
presence of this feature in members of the Themis family is well established, 
it is absent among the Lixiaohua family members. 

Despite the differences with respect to the Themis family and lack of iden- 
tification of aqueous alteration traces in spectra of its members, two MBCs 
(313P and 358P) are robustly associated with the Lixiaohua family. To better 
understand the family and its links to the MBCs, further analyses are needed. 
In particular, observations at the 3 ¡um region would be helpful to clarify if 
hydrated minerals are indeed not present in family members, as suggested 
based on the 0.7 ¡um band. 


Theobalda family.- This is the youngest and smallest of three families with 
two or more MBCs associated. The Theobalda family was formed by an impact 
cratering on the largest fragment, asteroid (778) Theobalda, that is found to 
occur about 7+ 2 Myr ago (Novaković, 2010). The author also found that the 
family originates from a parent body that was about 80-100 km in size. This 
is a significantly smaller parent body than in the Themis family's case and 
probably about factor 2 smaller than the parent body of the Lixiaohua family. 

The largest member of the family, (778) Theobalda, has been spectroscopi- 
cally classified as an F-type asteroid (Tholen, 1989). The data on other family 
members are limited, but available taxonomic classifications include C-, F-, 
and X-type objects, indicating their likely primitive compositions. This is fur- 
ther supported with an average albedo of family members of about p, — 0.062, 
based on the data of about 100 objects (Masiero et al., 2015). Apart from 
this limited data, Theobalda is a relatively poorly studied family, and addi- 
tional observations are strongly encouraged to characterise the family members 
better. 

Two MBCs (P/2016 J1 and 427P) are dynamically associated with the 
family. For both objects, sublimation of water-ice seems to be the most plausi- 
ble mechanism of the mass loss. Nevertheless, both MBCs are very intriguing 
and deserve additional studies. The P/2016 J1 is actually a double component 
object, and in this respect, we can also say that there are three MBCs asso- 
ciated with the Theobalda family. Moreno et al. (2017) found that P/2016 J1 
split into two components (J1-A and J1-B) only several years ago, making it 
also the youngest known the so-called asteroid pair. Yet, the authors found 
that the separation event and the present dust activity are unrelated. The 
427P (P/2017 S5) is interesting from another point of view. Its time-series pho- 
tometry lightcurve indicates possibly rapid rotation with a period of 1.4 hr. 
Though this period is yet to be confirmed, Jewitt et al. (2019) suggested that, 
if real, it may also play a role in the mass loss. The authors further suggested 
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Fig. 13 Members of three outer main-belt asteroid families containing at least two main- 
belt comets projected into the space of proper orbital elements. The locations of associated 
MBGs are also indicated with red circles. The background asteroids are shown as gray dots. 


that such a rapid rotation could be a consequence of spin-up by sublimation 
torques. These findings open new questions on the relation and relative impor- 
tance of sublimation and rotationally induced mass loss. For instance, one 
possibility is that the fast rotation is preceding the sublimation-driven activity, 
possibly causing surface landslides that may expose at the surface some shal- 
low ice, which in turn triggers the sublimation driven activity. Alternatively, 
the sublimation-driven activity produces torques that could spin up the small 
asteroids to the point when they also start to shed some mass due to the fast 
rotation or even disintegrate into smaller fragments. 
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4.3.2 Future prospects and potentially new reservoirs of 
main-belt comets 


Interestingly, the three families with the largest number of associated MBCs 
are of quite different ages. Even considering that 133P and 288P could be 
members of much younger subfamilies of the Themis family, the young age 
does not seem to be a key for a family to contain MBCs. Therefore, it seems 
that links between the families and main-belt comets are primarily based on 
family compositions and not on their ages. In this respect, the dark families of 
primitive composition are the main reservoirs of MBCs, and, therefore, could 
also represent large water-ice reservoirs in the main asteroid belt. 

There is still a lot of work to be done to understand these links prop- 
erly. De Prá et al. (2020a) recently studied the diversity of physical properties 
among five primitive outer main-belt families: Themis, Hygiea, Ursula, Veri- 
tas, and Lixiaohua families. The authors found that members of the Themis 
and Hygiea families show similar hydration levels. In contrast, the Ursula and 
Lixiaohua families present no sign of hydrated members based on the analysis 
of visible spectra. The Veritas family, however, presents the highest fraction of 
hydrated members. These results are difficult to interpret based on the current 
associations of main-belt comets and families. In particular, why the Lixiaohua 
family, which is undoubtedly linked to some MBCs, shows no signs of hydra- 
tion, while members of the Veritas family, still not linked to any MBCs, seem 
to be the most hydrated objects among the five studied primitive families? 
Certainly, it could be because many MBCs are yet to be discovered, and we 
might find many of them among the Veritas family members. A good example 
is the first association of an active asteroid, namely P/2017 S8, to the Pallas 
family (Novakovié and Hsieh, 2019). 

A potentially significant reservoir of MBCs could also be the Euphrosyne 
asteroid family, yet another group of primitive C-type asteroids in the outer 
main belt (Yang et al., 2020b). The recent work by Yang et al. (2020a) sug- 
gested that the largest family member, asteroid (31) Euphrosyne, probably 
contains a significant fraction of water ice in its interior. This further suggests 
that the other family members could contain some water-ice. There is, however, 
still no known main-belt comet belonging to this family, though P/2016 P1 
(PANSTARRS) could be a member of the Euphrosyne family according to 
Novakovié and Hsieh (2019). 

Many new MBCs have been discovered every year, and this trend could 
even accelerate in the near future. This will open numerous opportunities to 
study further links between the families and main-belt comets. 

Let us conclude this section with another interesting point of some MBCs 
being potentially binary asteroids. Agarwal et al. (2017) found that the active 
member of the 288P group is a binary asteroid. Additionally, one of the 
MBCs belonging to the Theobalda family (see below), namely P/2016 J1 
(PANSTARRS), is a recently formed asteroid pair, which could also be previ- 
ously a binary object. More recently, Jewitt et al. (2021) obtained a lightcurve 
of another MBC, namely 331P /Gibbs (P/2012F5), and found indications that 
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it could be a contact binary asteroid. These results open questions about 
whether there is a link between the MBCs and binary asteroid formation and 
how both these populations are related to asteroid families. Note that a binary 
MBC belongs to a primitive young asteroid family in each of these three cases. 
On the other hand, some members of the very young Datura family could 
be contact binaries as well (see Section 4.1), but the Datura family is not of 
a primitive composition. While we do expect contact binaries to form in the 
process of re-accumulation of fragments following asteroid disruption (Campo 
Bagatin et al., 2020), it remains to be explained how the MBCs fit into the 
more global picture. 


5 Summary and Conclusions 


Asteroid families-related topics are an emerging field of solar system research. 
More than a century of investigation of these ice-rocky relatives brought a 
lot of big discoveries and outstanding accomplishments. All of these help us 
better understand the families and the solar system as a whole. Nevertheless, 
there are still open questions and problems, and the future is bringing new 
challenges. Ongoing and future ground- and space-based instruments supply 
new observational data that is increasing fast in qualitative and quantitive 
ways. The large data sets will raise additional challenges for scientists working 
on asteroid families. Likely, big challenges always create big opportunities, and 
we anticipate a bright future for asteroid families related research. 

A large number-density of asteroids will require substitute techniques to 
identify new asteroid families or attribute recent members to the existing 
ones. Some promising approaches are already in place. Nevertheless, these 
techniques will need further development and improvement to cope with the 
rapid increase of new asteroid discoveries. Combining different methods such 
as multi-step and machine learning-based algorithms could be one promising 
direction to follow. In this respect, it is worth reminding that any development 
involving machine learning should also include steps necessary to understand 
the algorithms and secure the reproducibility of the results. 

Age determination techniques improved significantly in the last about 15 
years, primarily through the backward orbit integration methods for young 
families and Yarkovsky-based chronology for older ones. However, the uncer- 
tainties in their determination, especially for old families, are still often 
considerable, and therefore this issue needs to be addressed soon. The first step 
in this direction could be extensive testing of existing methods. On one side, 
this should allow quantifying their uncertainties better. On the other hand, 
testing should identify the principal sources of errors, setting foundations for 
further improvements in age determination methods. 

Studies of older families (>50 Myr old) will benefit more from the new 
physical data and refined evolution models. Combined with more accurate 
ages, these will allow a more authentic and detailed reconstruction of the 
long-term evolution. The area of young asteroid families («10 Myr old), and 
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especially very young ones («1.5 Myr), is already advancing with the discovery 
of smaller asteroids. Many of these groups are increasing fast in the number 
of known members, allowing much deeper analysis. This will be the case even 
more in the near future. Additionally, the number of young families is far 
from complete, with many of them waiting to be discovered. This will increase 
the statistical sample of young families and provide a view of the diversity of 
the population. We recently already saw an example of such progress when 
Vokrouhlicky et al. (2021b) found that the Hobson family could be formed 
from the binary parent body. Nevertheless, this is just the tip of the iceberg. 
Many discoveries are yet to happen. 

Ultimately, revealing the secrets of asteroid families is not a single-step 
process but rather an accumulation of knowledge over the years. Even the best 
and the most remarkable results represent, in principle, only the best we can 
learn about the families and their asteroids at the moment. The new data and 
better models will certainly, in the future, push our understanding of families 
even beyond our current imaginations. 
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